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We seek to derive the structure of hadrons from the fundamental theory of the 



strong interaction, QCD. Our work is founded on the hypothesis that a constituent 
approximation can be derived from QCD, so that a relatively small number of quark 
and gluon degrees of freedom need be explicitly included in the state vectors for low- 
lying hadrons. 1 To obtain a constituent picture, we use a Hamiltonian approach in 
light-front coordinates. 2 

I do not believe that light-front Hamiltonian field theory is extremely useful for the 
study of low energy QCD unless a constituent approximation can be made, and I do not 
Ph. believe such an approximation is possible unless cutoffs that it violate manifest gauge 

invariance and covariance are employed. Such cutoffs inevitably lead to relevant and 
marginal effective interactions {i.e., counterterms) that contain functions of longitudinal 
momenta. 3, 4 It is it not possible to renormalize light-front Hamiltonians in any useful 
manner without developing a renormalization procedure that can produce these non- 
canonical counterterms. 

The line of investigation I discuss has been developed by a small group of theorists 
who are working or have worked at Ohio State University and Warsaw University. 1 ' 3 " 19 
Ken Wilson provided the initial impetus for this work, and at a very early stage outlined 
much of the basic strategy we employ. 15 

I make no attempt to provide enough details to allow the reader to start doing 
light-front calculations. The introductory article by Harindranath 16 is helpful in this 
regard. An earlier version of these lectures 8 also provides many more details. 

A Constituent Approximation Depends on Tailored Renormalization 

If it is possible to derive a constituent approximation from QCD, we can formulate 
the hadronic bound state problem as a set of coupled few-body problems. We obtain 
the states and eigenvalues by solving 

H A \V A ) = E\V A ), (1) 

where, 

I * a) = < 5 I W) + I M9) + ■■■■ (2) 
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Figure 1. The second-order shift in the self-energy of a bare electron due to mixing with 
electron-photon states. 

where I use shorthand notation for the Fock space components of the state. The full 
state vector includes an infinite number of components, and in a constituent approx- 
imation we truncate this series. We derive the Hamiltonian from QCD, so we must 
allow for the possibility of constituent gluons. I have indicated that the Hamiltonian 
and the state both depend on a cutoff, A, which is critical for the approximation. 

This approach has no chance of working without a renormalization scheme tai- 
lored to it. Much of our work has focused on the development of such a renormalization 
scheme. 3 ' 4,17-19 In order to understand the constraints that have driven this develop- 
ment, seriously consider under what conditions it might be possible to truncate the 
above series without making an arbitrarily large error in the eigenvalue. I focus on 
the eigenvalue, because it will certainly not be possible to approximate all observable 
properties of hadrons (e.g., wee parton structure functions) this way. 

For this approximation to be valid, all many-body states must approximately de- 
couple from the dominant few-body components. We know that even in perturbation 
theory, high energy many-body states do not simply decouple from few-body states. In 
fact, the errors from simply discarding high energy states are infinite. In second-order 
perturbation theory, for example, high energy photons contribute an arbitrarily large 
shift to the mass of an electron. This second-order effect is illustrated in Figure 1, and 
the precise interpretation for this light-front time-ordered diagram will be given below. 
The solution to this problem is well-known, renormalization. We must use renormaliza- 
tion to move the effects of high energy components in the state to effective interactions^] 
in the Hamiltonian. 

It is difficult to see how a constituent approximation can emerge using any regu- 
larization scheme that does not employ a cutoff that either removes degrees of freedom 
or removes direct couplings between degrees of freedom. A Pauli-Villars "cutoff," for 
example, drastically increases the size of Fock space and destroys the hermiticity of the 
Hamiltonian. 

In the best case scenario we expect the cutoff to act like a resolution. If the cutoff 
is increased to an arbitrarily large value, the resolution increases and instead of seeing 
a few constituents we resolve the substructure of the constituents and the few-body 
approximation breaks down. As the cutoff is lowered, this substructure is removed 
from the state vectors, and the renormalization procedure replaces it with effective 

* These include one-body operators that modify the free dispersion relations. 
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interactions in the Hamiltonian. Any "cutoff" that does not remove this substructure 
from the states is of no use to us. 

This point is well-illustrated by the QED calculations discussed below 8 ' 10 ' n . There 
is a window into which the cutoff must be lowered for the constituent approximation to 
work. If the cutoff is too large, atomic states must explicitly include photons. After the 
cutoff is lowered to a value that can be self-consistently determined a-posteriori, photons 
are removed from the states and replaced by the Coulomb interaction and relativistic 
corrections. The cutoff cannot be lowered too far using a perturbative renormalization 
group, hence the window. 

Thus, if we remove high energy degrees of freedom, or coupling to high energy 
degrees of freedom, we should encounter self-energy shifts leading to effective one-body 
operators, vertex corrections leading to effective vertices, and exchange effects lead- 
ing to explicit many-body interactions not found in the canonical Hamiltonian. We 
naively expect these operators to be local when acting on low energy states, because 
simple uncertainty principle arguments indicate that high energy virtual particles can- 
not propagate very far. Unfortunately this expectation is indeed naive, and at best we 
can hope to maintain transverse locality. 4 I will elaborate on this point below. The 
study of perturbation theory with the cutoffs we must employ makes it clear that it 
is not enough to adjust the canonical couplings and masses to renormalize the theory. 
It is not possible to make significant progress towards solving light-front QCD without 
fully appreciating this point. 

Low energy many-body states do not typically decouple from low energy few- 
body states. The worst of these low energy many-body states is the vacuum. This 
is what drives us to use light-front coordinates. 2 Figure 2 shows a pair of particles 
being produced out of the vacuum in equal-time coordinates t and z. The transverse 
components x and y are not shown, because they are the same in equal-time and 
light-front coordinates. The figure also shows light-front time, 

x + = t + z , 

and the light-front longitudinal spatial coordinate, 

x~ = t - z . (4) 

In equal-time coordinates it is kinematically possible for virtual pairs to be pro- 
duced from the vacuum (although relevant interactions actually produce three or more 
particles from the vacuum), as long as their momenta sum to zero so that three- 
momentum is conserved. Because of this, the state vector for a proton includes an 
arbitrarily large number of particles that are disconnected from the proton. The only 
constraint imposed by relativity is that particle velocities be less than or equal to that 
of light. 

In light-front coordinates, however, we see that all allowed trajectories lie in the 
first quadrant. In other words, light-front longitudinal momentum, p + (conjugate to 
x~ since a ■ b = \(a + b~ + aTb + ) — aj_ ■ bj_), is always positive, 

P + > . (5) 

We exclude p + = 0, forcing the vacuum to be trivial because it is the only state with 
p + = 0. Moreover, the light-front energy of a free particle of mass m is 

p 2 ± +m 2 

P = t 

p+ 

3 



(3) 
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Figure 2. Light-front coordinates. Light-front 'time' is x + = t + z, and the light-front longitudinal 
spatial coordinate is x~ = t — z. 

This implies that all free particles with zero longitudinal momentum have infinite en- 
ergy, unless their mass and transverse momentum are identically zero. Replacing such 
particles with effective interactions should be reasonable. 

• Is the vacuum really trivial? 

• What about confinement? 

• What about chiral symmetry breaking? 

• What about instantons? 

• What about the job security of theorists who study the vacuum? 

The question of how one should treat "zero modes," degrees of freedom (which may 
be constrained) with identically zero longitudinal momentum, divides the light-front 
community. Our attitude is that explicitly including zero modes defeats the purpose 
of using light-front coordinates, and we do not believe that significant progress will be 
made in this direction, at least not in 3 + 1 dimensions. 

The vacuum in our formalism is trivial. We are forced to work in the "hidden 
symmetry phase" of the theory, and to introduce effective interactions that reproduce 
all effects associated with the vacuum in other formalisms. 1 ' 20 ' 21 The simplest exam- 
ple of this approach is provided by a scalar field theory with spontaneous symmetry 
breaking. It is possible to shift the scalar field and deal explicitly with a theory con- 
taining symmetry breaking interactions. In the simplest case 3 is the only relevant 
or marginal symmetry breaking interaction, and one can simply tune this coupling to 
the value corresponding to spontaneous rather than explicit symmetry breaking. Ken 
Wilson and I have also shown that in such simple cases one can use coupling coherence 
to fix the strength of this interaction so that tuning is not required. 3 
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I will make an additional drastic assumption in these lectures, an assumption that 
Ken Wilson does not believe will hold true. I will assume that all effective interactions 
we require are local in the transverse direction. If this is true, there are a finite number 
of relevant and marginal operators, although each contains a function of longitudinal 
momenta that must be determined by the renormalization procedure.^ There are many 
more relevant and marginal operators in the renormalized light-front Hamiltonian than 
in £qcd- If transverse locality is violated, the situation is much worse than this. 

The presence of extra relevant and marginal operators that contain functions 
tremendously complicates the renormalization problem, and a common reaction to 
this problem is denial, which may persist for years. However, this situation may make 
possible tremendous simplifications in the final nonperturbative problem. For example, 
few-body operators must produce confinement manifestly! 

Confinement cannot require particle creation and annihilation, flux tubes, etc. This 
is easily seen using a variational argument. Consider a color neutral quark- ant iquark 
pair that are separated by a distance R, which is slowly increased to infinity. Moreover, 
to see the simplest form of confinement assume that there are no light quarks, so that 
the energy should increase indefinitely as they are separated if the theory possesses 
confinement. At each separation the gluon components of the state adjust themselves 
to minimize the energy. But this means that the expectation value of the Hamiltonian 
for a state with no gluons must exceed the energy of the state with gluons, and there- 
fore must diverge even more rapidly than the energy of the true ground state. This 
means that there must be a two-body confining interaction in the Hamiltonian. If the 
renormalization procedure is unable to produce such confining two-body interactions, 
the constituent picture will not arise. 

Manifest gauge invariance and manifest rotational invariance require all physical 
states to contain an arbitrarily large number of constituents. Gauge invariance is not 
manifest since we work in light-cone gauge with the zero modes removed; and it is easy 
to see that manifest rotational invariance requires an infinite number of constituents. 
Rotations about transverse axes are generated by dynamic operators in interacting 
light-front field theories, operators that create and annihilate particles. 2, 22 No state 
with a finite number of constituents rotates into itself or transforms as a simple tensor 
under the action of such generators. These symmetries seem to imply that we can not 
obtain a constituent approximation. 

To cut this Gordian knot we employ cutoffs that violate gauge invariance and 
covariance, symmetries which then must be restored by effective interactions, and which 
need not be restored exactly. A familiar example of this approach is supplied by lattice 
gauge theory, where rotational invariance is violated by the lattice. 

Simple Strategy 

We have recently employed a conceptually simple strategy to complete bound 
state calculations. The first step is to use a perturbative similarity renormalization 
group 17 ' 18 ' 19 and coupling coherence 3 ' 4 to find the renormalized Hamiltonian as an 
expansion in powers of the canonical coupling: 

H A = h + g A h J t + glh% + --- (7) 

We compute this series to a finite order, and to date have not required any ad hoc as- 
sumptions to uniquely fix the Hamiltonian. No operators are added to the hamiltonian, 
so the hamiltonian is completely determined by the underlying theory to this order. 

^These functions imply that there are effectively an infinite number of relevant and marginal operators; 
however, their dependence on fields and transverse momenta is extremely limited. 
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The second step is to employ bound state perturbation theory to solve the eigen- 
value problem. The complete Hamiltonian contains every interaction (although each 
is cut off) contained in the canonical Hamiltonian, and many more. We separate the 
Hamiltonian, 

H A = H A + V A , (8) 

treating H A nonperturbatively and computing the effects of V A in bound state pertur- 
bation theory. We must choose H A and A so that H A is manageable and to minimize 
corrections from higher orders of V A within a constituent approximation. 

If a constituent approximation is valid after A is lowered to a critical value which 
must be determined, we may be able to move all creation and annihilation operators 
to V A . will include many-body interactions that do not change particle number, 
and these interactions should be primarily responsible for the constituent bound state 
structure. 

There are several obvious flaws in this strategy. Chiral symmetry-breaking op- 
erators, which must be included in the Hamiltonian since we work entirely in the 
hidden symmetry phase of the theory, do not appear at any finite order in the cou- 
pling. These operators must simply be added and tuned to fit spectra or fixed by a 
non-perturbative renormalization procedure. 1 ' 7 ' 9 In addition, there are perturbative 
errors in the strengths of all operators that do appear. We know from simple scaling 
arguments 23 that when A is in the scaling regime: 

• small errors in relevant operators exponentiate in the output, 

• small errors in marginal operators produce comparable errors in output, 

• small errors in irrelevant operators tend to decrease exponentially in the output. 

This means that even if a relevant operator appears {e.g., a constituent quark or gluon 
mass operator), we may need to tune its strength to obtain reasonable results. We have 
not had to do this, but we have recently studied some of the effects of tuning a gluon 
mass operator. 14 

To date this strategy has produced well-known results in QED 8, 10 ' 11 through the 
Lamb shift, and reasonable results for heavy quark bound states in QCD. 8 ' 12, 13 ' 14 The 
primary objective of the remainder of these lectures is to review these results. I first 
use the Schwinger model as an illustration of a light-front bound state calculation. This 
model does not require renormalization, so before turning to QED and QCD I discuss 
the renormalization procedure that we have developed. 



LIGHT-FRONT SCHWINGER MODEL 

In this section I use the Schwinger model, massless QED in 1 + 1 dimensions, 
to illustrate the basic strategy we employ after we have computed the renormalized 
Hamiltonian. No model in 1 + 1 dimensions illustrates the renormalization problems 
we must solve before we can start to study QCD 3+1 . 

The Schwinger model can be solved analytically. 24 ' 25 Charged particles are con- 
fined because the Coulomb interaction is linear and there is only one physical particle, 
a massive neutral scalar particle with no self-interactions. The Fock space content of 
the physical states depends crucially on the coordinate system and gauge, and it is only 
in light-front coordinates that a simple constituent picture emerges. 26 
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The Schwinger model was first studied in Hamiltonian light-front field theory by 
Bergknoff. 26 My description of the model follows his closely, and I recommend his pa- 
per to the reader. Bergknoff showed that the physical boson in the light-front massless 
Schwinger model in light-cone gauge is a pure electron-positron state. This is an amaz- 
ing result in a strong-coupling theory of massless bare particles, and it illustrates how 
a constituent picture may arise in QCD. The electron-positron pair is confined by the 
linear Coulomb potential. The light-front kinetic energy vanishes in the massless limit, 
and the potential energy is minimized by a wave function that is flat in momentum 
space, as one might expect since a linear potential produces a state that is as localized 
as possible (given kinematic constraints due to the finite velocity of light) in position 
space. 

In order to solve this theory I must first set up a large number of details. I 
recommend that for a first reading these details be skimmed, because the general idea 
is more important than the detailed manipulations. The Lagrangian for the theory is 

C = $(i?- m)^j - ^ lv $A* - \ F ^ V , (9) 

where F^ u is the electromagnetic field strength tensor. I have included an electron 
mass, m, which is taken to zero later. I choose light-cone gauge, 

A + = . (10) 

In this gauge we avoid ghosts, so that the Fock space has a positive norm. This is 
absolutely essential if we want to apply intuitive techniques from many-body quantum 
mechanics. 

Many calculations are simplified by the use of a chiral representation of the Dirac 
gamma matrices, so in this section I will use: 





(11) 

which leads to the light-front coordinate gamma matrices, 



In light-front coordinates the fermion field ip contains only one dynamical degree 
of freedom, rather than two. To see this, first define the projection operators, 

^°Hoo) a-sMs:)- (13) 

Using these operators split the fermion field into two components, 

i); = + ^_ = + A_?/> • (14) 
The two-component Dirac equation in this gauge is 

{b +d ~ + b~ d+ ~ m ~ h +A ~)^ = ; (15) 

which can be split into two one-component equations, 

id~ip + = —imip- + eA~ip + , (16) 
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id + ip- = irmp + . (17) 

Here if}± refers to the non-zero component of ip±. 

The equation for ifj + involves the light-front time derivative, d~; so ip + is a dy- 
namical degree of freedom that must be quantized. On the other hand, the equation 
for ip_ involves only spatial derivatives, so is a constrained degree of freedom that 
should be eliminated in favor of Formally, 

1>- = ^+ • (18) 

This equation is not well-defined until boundary conditions are specified so that d + can 
be inverted. I will eventually define this operator in momentum space using a cutoff, 
but I want to delay the introduction of a cutoff until a calculation requires it. 
I have chosen the gauge so that A + = 0, and the equation for A~ is 

- \(d + ?A- = e*pU + . (19) 

A" is also a constrained degree of freedom, and we can formally eliminate it, 

4e 



(«9+) 



• (20) 



We are now left with a single dynamical degree of freedom, ip + , which we can 
quantize at x + = 0, 

{Mx-),^ + (y-)} = A + 8(x--y-). (21) 

We can introduce free particle creation and annihilation operators and expand the field 
operator at x + = 0, 



dk+ 



lk+>o 4tt 
with, 



b k e~ ik - x + d{e ik - x 



(22) 



{b k ,bl} = 4n5(k+-p+). (23) 

In order to simplify notation, I will often write k to mean k + . If I need k~ = m 2 /k + , I 
will provide the superscript. 

The next step is to formally specify the Hamiltonian. I start with the canonical 
Hamiltonian, 

H = j dx~(H + v) , (24) 



id+ 



(26) 



To actually calculate: 

• replace ip + with its expansion in terms of b k and d k , 

• normal-order, 

• throw away constants, 



8 



• drop all operators that require 60 and do- 
The free part of the Hamiltonian becomes 



When V is normal-ordered, we encounter new one-body operators, 



e 2 r dk 
2ir Jk>o 47r 



p>o^ P ((k - p) 2 



1 



ik+p)' 



(b\b k + 44) • 



(27) 



(2? 



This operator contains a divergent momentum integral. From a mathematical point of 
view we have been sloppy and need to carefully add boundary conditions and define 
how d + is inverted. However, I want to apply physical intuition and even though no 
physical photon has been exchanged to produce the initial interaction, I will act as if a 
photon has been exchanged and everywhere an 'instantaneous photon exchange' occurs 
I will cut off the momentum. In the above integral I insist, 



\k — p\ > e . 



Using this cutoff we find that 

e 2 rdk/1 
it J Air\e 



( l --^ + 0(e))(b\b k + dld k ). 



1 

k 



(29) 



(30) 



Comparing this result with the original free Hamiltonian, we see that a divergent 
mass-like term appears; but it does not have the same dispersion relation as the bare 
mass. Instead of depending on the inverse momentum of the fermion, it depends on the 
inverse momentum cutoff, which cannot appear in any physical result. There is also a 
finite shift in the bare mass, with the standard dispersion relation. 

The normal-ordered interactions are 

dki dk 4 



V = 2e 2 f ^ • • • ^AnSih + k 2 - k 3 - k 4 ) 
J Att Ait 



b\d\d 4 b 3 



(h - fc 3 ) 



b\d\d 3 b 4 



(h + k 2 ) 



(h - k 3 y 



(b\blb 3 b 4 + d\d\d 3 d 4 ) + 



(31) 



I do not display the interactions that involve the creation or annihilation of elec- 
tron/positron pairs, which are important for the study of multiple boson eigenstates. 

The first term in V is the electron-positron interaction. The longitudinal momen- 
tum cutoff I introduced above requires \k\ — k 3 \ > e, so in position space we encounter 
a potential which I will naively define with a Fourier transform that ignores the fact 
that the momentum transfer cannot exceed the momentum of the state, 



v(x ) 



Q1Q2 



^ ^ 0(\k\ - e) expi-'-kx ) 



7T 



7T 

- —\x 



(9(e) 



(32) 



This potential contains a linear Coulomb potential that we expect in two dimensions, 
but it also contains a divergent constant that is negative for unlike charges and positive 
for like charges. 
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In charge neutral states the infinite constant in V is exactly canceled by the di- 
vergent 'mass' term in Hq. This Hamiltonian assigns an infinite energy to states with 
net charge, and a finite energy as e — > to charge zero states. This does not imply 
that charged particles are confined, but the linear potential prevents charged particles 
from moving to arbitrarily large separation except as charge neutral states. The con- 
finement mechanism I propose for QCD in 3+1 dimensions shares many features with 
this interaction. 

I would also like to mention that even though the interaction between charges 
is long-ranged, there are no van der Waals forces in 1+1 dimensions. It is a simple 
geometrical calculation to show that all long range forces between two neutral states 
cancel exactly. This does not happen in higher dimensions, and if we use long-range 
two-body operators to implement confinement we must also find many-body operators 
that cancel the strong long-range van der Waals interactions. 

Given the complete Hamiltonian in normal-ordered form we can study bound 
states. A powerful tool for the initial study of bound states is the variational wave 
function. In this case, we can begin with a state that contains a single electron-positron 
pair, 

rP dp 



MP)) = I £<t>{p)b ] / P _ p \V) . (33) 



The norm of this state is 



^(P')\^(P))=AnP5(P'-P)^J o P J|0(P)| 2 } , (34) 

where the factors outside the brackets provide a covariant plane wave normalization for 
the center-of-mass motion of the bound state, and the bracketed term should be set to 
one. 

The expectation value of the one-body operators in the Hamiltonian is 



{9\Ho + H' \*) = jf 



dk 
An 



m 2 — e 2 lix m 2 — e 2 /ir 2e 2 

+ ^ 7— + 



k P — k Tie 

and the expectation value of the normal-ordered interactions is 

4e 2 fdk 1 dk 2 r 1 1 



(35) 



4e r' 



47T 47T 



(h - k 2 ) 2 P 2 



rihMkz) , (36) 



where I have dropped the overall plane wave norm. The prime on the last integral 
indicates that the range of integration in which \ki — k 2 \ < e must be removed. By 
expanding the integrand about ki — k 2 , one can easily confirm that the 1/e divergences 
cancel. 

With m = the energy is minimized when 

(f){k) = , (37) 

and the invariant-mass is 

e 2 

M 2 = — . (38) 

7T 

This type of simple analysis can be used to show that this electron-positron state 
is actually the exact ground state of the theory with momentum P, and that bound 
states do not interact with one another. 

The primary purpose of introducing the Schwinger model is to illustrate that 
bound state center-of-mass motion is easily separated from relative motion in light- 
front coordinates, and that standard quantum mechanical techniques can be used to 
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analyze the relative motion of charged particles once the Hamiltonian is found. It is 
intriguing that even when the fermion is massless, the states are constituent states 
in light-cone gauge and in light-front coordinates. This is not true in other gauges 
and coordinate systems. The success of light-front field theory in 1+1 dimensions can 
certainly be downplayed, but it should be emphasized that no other method on the 
market is as powerful for bound state problems in 1+1 dimensions. 

The most significant barriers to using light-front field theory to solve low en- 
ergy QCD are not encountered in 1+1 dimensions. The Schwinger model is super- 
renormalizable, so we completely avoid serious ultraviolet divergences. There are no 
transverse directions, and we are not forced to introduce a cutoff that violates rotational 
invariance, because there are no rotations. Confinement results from the Coulomb in- 
teraction, and chiral symmetry is not spontaneously broken. This simplicity disappears 
in realistic 3 + 1-dimensional calculations, which is one reason there are so few 3 + 1- 
dimensional light-front field theory calculations. 



LIGHT-FRONT RENORMALIZATION GROUP: SIMILARITY TRANS- 
FORMATION AND COUPLING COHERENCE 

As argued above, in 3 + 1 dimensions we must introduce a cutoff on energies, 
A, and we never perform explicit bound state calculations with A anywhere near its 
continuum limit. In fact, we want to let A become as small as possible. In my opinion, 
any strategy for solving light-front QCD that requires the cutoff to explicitly approach 
infinity in the nonperturbative part of the calculation is useless. Therefore, we must 
set up and solve 

p 2 _i_ /if 2 

Pa I *a(P)> = p + I *a(P)> • (39) 

Physical results, such as the mass, M, can not depend on the arbitrary cutoff, 
A, even as A approaches the scale of interest. This means that P^ and | ^a) must 
depend on the cutoff in such a way that (^a I Pa I ^a) does n °t- Wilson based 
the derivation of his renormalization group on this observation, 27 ' 28 ' 29 ' 23 and we use 
Wilson's renormalization group to compute P^. 

It is difficult to even talk about how the Hamiltonian depends on the cutoff without 
having a means of changing the cutoff. If we can change the cutoff, we can explicitly 
watch the Hamiltonian's cutoff dependence change and fix its cutoff dependence by 
insisting that this change satisfy certain requirements (e.g., that the limit in which the 
cutoff is taken to infinity exists). We introduce an operator that changes the cutoff, 

#(Ai) = T[H(A )} , (40) 

where I assume that A 1 < A . To simplify the notation, I will let H(A[) = Hi. To 
renormalize the hamiltonian we study the properties of the transformation. 

Figure 3 displays two generic cutoffs that might be used. Traditionally theorists 
have used cutoffs that remove high energy states, as shown in Figure 3a. This is the type 
of cutoff Wilson employed in his initial work 27 and I have studied its use in light-front 
field theory. 4 When a cutoff on energies is reduced, all effects of couplings eliminated 
must be moved to effective operators. As we will see explicitly below, when these effec- 
tive operators are computed perturbatively they involve products of matrix elements 
divided by energy denominators. Expressions closely resemble those encountered in 
standard perturbation theory, with the second-order operator involving terms of the 
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Figure 3. Two ways to run a cutoff on free energy. In (a) a cutoff on the magnitude of the energy is 
lowered from the solid to the dashed lines, with problems resulting from the removed shaded region. 
In (b) a cutoff on how far off diagonal matrix elements appear is lowered from the dashed to the solid 
lines. 

form 

5Vij „ i v_ i mi i y i m . (4i) 

£i — Cfc 

This new effective interaction replaces missing couplings, so the states fa and <pj are 
retained and the state (f>k is one of the states removed. The problem comes from the 
shaded, lower right-hand corner of the matrix, where the energy denominator van- 
ishes for states at the corner of the remaining matrix. In this corner we should use 
nearly degenerate perturbation theory rather than perturbation theory, but to do this 
requires solving high energy many-body problems nonperturbatively before solving the 
low energy few-body problems. 

An alternative cutoff, which does not actually remove any states and which can be 
run by a similarity transformational is shown in Figure 3b. This cutoff removes couplings 
between states whose free energy differs by more than the cutoff. The advantage of this 
cutoff is that the effective operators resulting from it contain energy denominators which 
are never smaller than the cutoff, so that a perturbative approximation for the effective 
Hamiltonian may work well. I discuss a conceptually simple similarity transformation 
that runs this cutoff below. 

Given a cutoff and a transformation that runs the cutoff, we can discuss how the 
Hamiltonian depends on the cutoff by studying how it changes with the cutoff. Our 
objective is to find a "renormalized" Hamiltonian, which should give the same results 
as an idealized Hamiltonian in which the cutoff is infinite and which displays all of 
the symmetries of the theory. To state this in simple terms, consider a sequence of 
Hamiltonians generated by repeated application of the transformation, 

H - H x - H 2 - • • • . (42) 

What we really want to do is fix the final value of A at a reasonable hadronic scale, 
and let Ao approach infinity. In other words, we seek a Hamiltonian that survives an 

■■'In deference to the original work I will call this a similarity transformation even though in all cases 
of interest to us it is a unitary transformation. 
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infinite number of transformations. In order to do this we need to understand what 
happens when the transformation is applied to a broad class of Hamiltonians. 

Perturbative renormalization group analyses typically begin with the identification 
of at least one fixed point, H*. A fixed point is defined to be any Hamiltonian that 
satisfies the condition 

H* = T[H*] . (43) 

For perturbative renormalization groups the search for such fixed points is relatively 
easy. If H* contains no interactions (i.e., no terms with a product of more than two field 
operators), it is called Gaussian. If H* has a massless eigenstate, it is called critical. If 
a Gaussian fixed point has no mass term, it is a critical Gaussian fixed point. If it has 
a mass term, this mass must typically be infinite, in which case it is a trivial Gaussian 
fixed point. In lattice QCD the trajectory of renormalized Hamiltonians stays 'near' a 
critical Gaussian fixed point until the lattice spacing becomes sufficiently large that a 
transition to strong-coupling behavior occurs. If H* contains only weak interactions, 
it is called near- Gaussian, and we may be able to use perturbation theory both to 
identify H* and to accurately approximate 'trajectories' of Hamiltonians near H*. Of 
course, once the trajectory leaves the region of H* , it is generally necessary to switch 
to a non-perturbative calculation of subsequent evolution. 

Consider the immediate 'neighborhood' of the fixed point, and assume that the 
trajectory remains in this neighborhood. This assumption must be justified a posteriori, 
but if it is true we should write 

H l = H* + 5H X , (44) 

and consider the trajectory of small deviations 5Hi. 

As long as 5Hi is 'sufficiently small,' we can use a perturbative expansion in powers 
of 5 Hi, which leads us to consider 

8H l+1 = L ■ 5Hi + NtfHi] . (45) 

Here L is the linear approximation of the full transformation in the neighborhood of 
the fixed point, and N[5Hi] contains all contributions to 5Hi + i of 0(5H?) and higher. 

The object of the renormalization group calculation is to compute trajectories and 
this requires a representation for 5 Hi. The problem of computing trajectories is one 
of the most common in physics, and a convenient basis for the representation of 5Hi 
is provided by the eigenoperators of L, since L dominates the transformation near the 
fixed point. These eigenoperators and their eigenvalues are found by solving 

L ■ O m = \ m O m . (46) 

If H* is Gaussian or near-Gaussian it is usually straightforward to find L, and its eigen- 
operators and eigenvalues. This is not typically true if H* contains strong interactions, 
but in QCD we hope to use a perturbative renormalization group in the regime of 
asymptotic freedom, and the QCD ultraviolet fixed point is apparently a critical Gaus- 
sian fixed point. For light-front field theory this linear transformation is a scaling of 
the transverse coordinate, the eigenoperators are products of field operators and trans- 
verse derivatives, and the eigenvalues are determined by the transverse dimension of 
the operator. All operators can include both powers and inverse powers of longitudinal 
derivatives because there is no longitudinal locality. 4 

Using the eigenoperators of L as a basis we can represent 5 Hi, 

SH t = Hmftm + 9m,O m + J2 W m l O m • (47) 

rneR rneM m£l 
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Here the operators O m with m G R are relevant (i.e., X m > 1), the operators O m with 
m G M are marginal (i.e., X m = 1), and the operators with m G / are either irrelevant 
(i.e., X m < 1) or become irrelevant after many applications of the transformation. The 
motivation behind this nomenclature is made clear by considering repeated application 
of L, which causes the relevant operators to grow exponentially, the marginal operators 
to remain unchanged in strength, and the irrelevant operators to decrease in magnitude 
exponentially. There are technical difficulties associated with the symmetry of L and 
the completeness of the eigenoperators that I ignore. 30 

For the purpose of illustration, let me assume that A m = 4 for all relevant operators, 
and A m = 1/4 for all irrelevant operators. The transformation can be represented by 
an infinite number of coupled, nonlinear difference equations: 

A*m; +1 = 4/i m; + [fjj mi , g mi , w mi ] , (48) 

9m l+1 = 9m l + Ng m [/i m; , g mv w rrH ] , (49) 

w mi+l = ^w mi + N Wm [fi mi , g mi , w mi \ . (50) 

Sufficiently near a critical Gaussian fixed point, the functions N^ m , N gm , and N Wm 
should be adequately approximated by an expansion in powers of fj, mt , g mi , and w mr 
The assumption that the Hamiltonian remains in the neighborhood of the fixed point, 
so that all \i mi , g mi , and w mi remain small, must be justified a posteriori. Any precise 
definition of the neighborhood of the fixed point within which all approximations are 
valid must also be provided a posteriori. 

Wilson has given a general discussion of how these equations can be solved, 23 but 
I will use coupling coherence 3 ' 4 to fix the Hamiltonian. This is detailed below, so at 
this point I will merely state that coupling coherence allows us to fix all couplings 
as functions of the canonical couplings and masses in a theory. The renormalization 
group equations specify how all of the couplings run, and coupling coherence uses this 
behavior to fix the strength of all of the couplings. But the first step is to develop a 
transformation. 



Similarity Transformation 

Stan Glazek and Ken Wilson studied the problem of small energy denominators 
which are apparent in Wilson's first complete non-perturbative renormalization group 
calculations, 28 and realized that a similarity transformation which runs a different form 
of cutoff (as discussed above) avoids this problem. 17 ' 18 ' 19 Independently, Wegner 32 
developed a similarity transformation which is easier to use than that of Glazek and 
Wilson. 

In this section I want to give a simplified discussion of the similarity transformation, 
using sharp cutoffs that must eventually be replaced with smooth cutoffs which require 
a more complicated formalism. 

Suppose we have a Hamiltonian, 

H A » = H A ° + V Ao , (51) 

where Hq° is diagonal. The cutoff A indicates that (4>i\V A °\(f)j} = if \E 0i — E j\ > A . 
I should note that A is defined differently in this section from later sections. We 
want to use a similarity transformation, which automatically leaves all eigenvalues and 
other physical matrix elements invariant, that lowers this cutoff to Ai. This similarity 
transformation will constitute the first step in a renormalization group transformation, 
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with the second step being a rescaling of energies that returns the cutoff to its original 
numerical value.| 

The transformed Hamiltonian is 



(H A ° + V A °)e 



(52) 



where R is a hermitian operator. If H A ° is already diagonal, R = 0. Thus, if R has an 
expansion in powers of V, it starts at first order and we can expand the exponents in 
powers of R to find the perturbative approximation of the transformation. 

We must adjust R so that the matrix elements of H Al vanish for all states that 
satisfy Ai < \E 0i — E j\ < A . We insist that this happens to each order in perturbation 
theory. Consider such a matrix element, 



\H 



Ai I 



z iR {H A ° + V A °) 



l) 



l + iR+-- -)(H A ° + V Ao )(l -iR )\(t>j) 



A > 



h\Ht\<l>i) + (<t>i\V A °+i 



R, Hq° 



/> + • 



(53) 



The last line contains all terms that appear in first-order perturbation theory. Since 
(<pi\H °\(/)j) = for these off-diagonal matrix elements, we can satisfy our new con- 
straint using 



Eoj — E ( 



Oi 



(54) 



E j\ < A to first order in V A °. 



E j I < Ai are zero to first order 



This fixes the matrix elements of R when Ai < \E 0i 
I will assume that the matrix elements of R for \E 0i 
in V, and fix these matrix elements to second order below. 

Given R we can compute the nonzero matrix elements of H Al . To second order in 
V A ° these are 



rrAi 
n ab 



^\H + V + i[R,Ho\ +i[R,V 
J a |ifo + V + i[i22,#ol|&> 



UR 2 ,H } + RH R + O(V 3 )\<p b ) 



- ^2 ®akQkbV a kVkb 
1 k 



-E 



e afc (i-e 



Eok — Eo a 

VakYkb 



+ 



1 



kb. 



Eok — Ei 



0a 



H + V + i R 2 ,H \cj> b ) 



+ - ®ak®kbVakVkb 
1 k 



Eok — E, 



Eok 

+ e fc6 (i - 
+ 



- E Qb 

®ak) 



VakYkb 



Eok — Ei 



1 



o« 



E() a — Enk\ — \E[ 



Ok 



>06 



Eok\) 



Eok — E ob 
-6(\E ob 



X 



06 



0{V 3 



Eok\ — \E( 



Oa 



Eok\) 



^ VakVkb 
k 







ak ' 



em 



Oa 



E 



Ok I 



\E 



Ob 



Eok\) 







Eok ~ E Qa 
\Eob — E Q k\ ■ 



\E< 



Oa 



kb- 



+ 

Eok\) 



E 



Ok 



E 



Ob 



(55) 



I have dropped the A superscript on the right-hand side of this equation and used 
subscripts to indicate matrix elements. The operator 0jj is one if Ai < \E Qi — E j\ < A 

§The rescaling step is not essential, but it avoids exponentials of the cutoff in the renormalization 
group equations. 
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and zero otherwise. It should also be noted that Vij is zero if IE* 



E 0j \ > A . All 

energy denominators involve energy differences that are at least as large as Ai, and this 
feature persists to higher orders in perturbation theory; which is the main motivation 
for choosing this transformation. R 2 is second-order in V, and we are still free to 
choose its matrix elements; however, we must be careful not to introduce small energy 
denominators when choosing R 2 . 

The matrix element {<j) a \i R2,Ho \4>b) = i(Eob — Eo a )((j) a \R 2 \<f) b ) must be specified. 
I will choose this matrix element to cancel the first sum in the final right-hand side of 
Eq. (55). This choice leads to the same result one obtains by integrating a differential 
transformation that runs a step function cutoff. To cancel the first sum in the final 
right-hand side of Eq. (55) requires 



{4>a\Rl\4>b) = T^^l^ak^kbVakVkb 



Eok — E, 



+ 



Oa 



E, 



Ok 



E 



Ob 



X 



0(\E Oa - E ok \ - \E ob - E ok \) - 9(\E ob - E ok \ - \E 0a - E { 



0k\ 



(56) 



No small energy denominator appears in R 2 because it is being used to cancel a term 
that involves a large energy difference. If we tried to use R 2 to cancel the remaining 
sum also, we would find that it includes matrix elements that diverge as E 0a — E ob goes 
to zero, and this is not allowed. 

The non-vanishing matrix elements of H Al are now completely determined to 



rrAi _ 
n ab — 



K\H + V 
^ VakVkb 



Ao I 



Q(\E 0a — E ok \ — \E ob — E ok \) 

Vak ^ ^ H 



e 



Eo k — E 0a 
0(\E ob — E ok \ — \E 0a — E ok \) 



Eo k — E, 



Ob 



(57) 



" A ° is zero if l^oi 



E j\ > Aq, so there are implicit cutoffs 



Let me again mention that 
that result from previous transformations. 

As a final word of caution, I should mention that the use of step functions produces 
long-range pathologies in the interactions that lead to infrared divergences in gauge 
theories. We must replace the step functions with smooth functions to avoid this 
problem. This problem will not show up in any calculations detailed in these lectures, 
but it does affect higher order calculations in QED 10 ' 11 and QCD. 31 



Light-Front Renormalization Group 



In this section I use the similarity transformation to form a perturbative light-front 
renormalization group for scalar field theory. If we want to stay as close as possible to 
the canonical construction of field theories, we can: 

• Write a set of 'allowed' operators using powers of derivatives and field operators. 

• Introduce 'free' particle creation and annihilation operators, and expand all field 
operators in this basis. 

• Introduce cutoffs on the Fock space transition operators. 
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Instead of following this program I will skip to the final step, and simply write a 
Hamiltonian to initiate the analysis. 



H = J dqi dq 2 (167r 3 )5 3 (gi - q 2 ) u 2 (-q u q 2 ) a\q 1 )a(q 2 ) 

+ ^ J dqi dq 2 dq 3 dq 4 (167r 3 )6 3 (<?i + q 2 + q 3 - q 4 )Q(q^ - q 3 - q 2 ~ <h) 

«4(-9i, -92, -Q3, q4> a\q 1 )a t (q 2 )a\q 3 )a(q i ) 

+- J dqi dq 2 dq 3 dq 4 (16vr 3 )5 3 (gi + q 2 - q 3 - 9 4 )6(gJ + q 3 ~ Q 2 ~ 9i~) 

Ui(-qi, -92, 93, 94) a ] {qi)a ] {q 2 )a(q 3 )a{q A ) 

+- J dq x dq 2 dq 3 dq 4 (167r 3 )5 3 (gi - q 2 - q 3 - 9 4 )6(g^ + q 3 + q 2 - 9i~) 

m 4 (-9i, 92, 93, 94) a t (gi)a(g 2 )a(93)a(94) 
+ 0{<f) , (58) 



where, 



dq + d 2 q L 

dq = Te^v ' (59) 



9L 



(60) 



and, 



e(Q-) = e(£-\Q-\) ■ (61) 



I assume that no operators break the discrete — > —0 symmetry. The functions 
u 2 and m 4 are not yet determined. If we assume locality in the transverse direction, 
these functions can be expanded in powers of their transverse momentum arguments. 
Note that to specify the cutoff both transverse and longitudinal momentum scales 
are required, and in this case the longitudinal momentum scale is independent of the 
particular state being studied. Note also that P + breaks longitudinal boost invariance 
and that a change in P + can be compensated by a change in A 2 . This may have 
important consequences, because the Hamiltonian should be a fixed point with respect 
to changes in the cutoff's longitudinal momentum scale, since this scale invariance is 
protected by Lorentz covariance. 4 

I have specified the similarity transformation in terms of matrix elements, and will 
work directly with matrix elements, which are easily computed in the free particle Fock 
space basis. In order to study the renormalization group transformation I will assume 
that the Hamiltonian includes only the interactions shown above. A single transforma- 
tion will produce a Hamiltonian containing products of arbitrarily many creation and 
annihilation operators, but it is not necessary to understand the transformation in full 
detail. 

I will define the full renormalization group transformation as: (i) a similarity 
transformation that lowers the cutoff in 0, (ii) a rescaling of all transverse momenta 
that returns the cutoff to its original numerical value, (iii) a rescaling of the creation 
and annihilation operators by a constant factor (, and (iv) an overall constant rescaling 
of the Hamiltonian to absorb a multiplicative factor that results from the fact that it 
has the dimension of transverse momentum squared. These rescaling operations are 
introduced so that it may be possible to find a fixed point Hamiltonian that contains 
interactions. 23 
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To find the critical Gaussian fixed point we need to study the linearized approx- 
imation of the full transformation, as discussed above. In general the linearized ap- 
proximation can be extremely complicated, but near a critical Gaussian fixed point 
it is particularly simple in light-front field theory with zero modes removed, because 
tadpoles are excluded. We have already seen that the similarity transformation does 
not produce any first order change in the Hamiltonian (see Eq. (57)), so the first order 
change is determined entirely by the rescaling operation. If we let 

Ai = V^o , (62) 

Pi± = VPi± , (63) 

and, 

a p = (a p , , a\ = (a^ , (64) 
to first order the transformed Hamiltonian is 



H = 



-12 ~Qi 



C 2 J dq x dq 2 (167T 3 )<5 3 (gi - q 2 ) u 2 {qf,mi) ^(q 1 )a(q 2 ) 

+ J dqi dq 2 dq 3 dq 4 (167r 3 )5 3 (gi + q 2 + q 3 - g 4 )0(<?4 - q 3 

Ui{qf, ml) a\q 1 )a\q 2 )a ] (q 3 )a(q 4 ) 

+ ^V 4 C 4 J dqi dq 2 dq 3 dq A (167r 3 )5 3 (gi + q 2 - q 3 - q^)<$>{ql + q 3 - q 2 - gf) 

Ui{qf, ml) a t (gi)a t (g 2 )a(g3)a(g 4 ) 

+ g^ 4 C 4 J dqi dq 2 dq 3 dq 4 (167r 3 )5 3 (gi - q 2 - q 3 - g 4 )©(g 4 ' + <h + <h ~ Qi) 

Ui{qf, mi) ^ (qi)a(q 2 )a(q 3 )a(q 4 ) 

+ W 6 ) , (65) 

where I have simplified my notation for the arguments appearing in the functions Ui. 
An overall factor of rj 2 that results from the engineering dimension of the Hamiltonian 
has been removed. The Gaussian fixed point is found by insisting that the first term 
remains constant, which requires 

( 2 u 2 (q + , i]q ± ) = u* 2 (q + , q ± ) . (66) 

I have used the fact that u 2 actually depends on only one momentum other than the 
cutoff. 

The solution to this equation is a monomial in which depends on (, 

/ -I \ (n/2) 

u* 2 (q + ,^) = f(q + )(^T, C=[-J • (67) 

The solution depends on our choice of n and to obtain the appropriate free particle 
dispersion relation we need to choose n = 2, so that 

^(<7 + ,q±) = /(<? + )ql , C = -- (68) 

V 

f(q + ) is allowed because the cutoff scale P + allows us to form the dimensionless variable 
q + 1 ' P + , which can enter the one-body operator. We will see that this happens in QED 
and QCD. Note that the constant in front of each four-point interaction becomes one, 
so that their scaling behavior is determined entirely by w 4 . If we insist on transverse 
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locality (which may be violated because we remove zero modes), we can expand u 4 
in powers of its transverse momentum arguments, and discover powers of rjq^ in the 
transformed Hamiltonian. Since we are lowering the cutoff, 77 < 1, and each power of 
transverse momentum will be suppressed by this factor. This means increasing powers 
of transverse momentum are increasingly irrelevant. 

I will not go through a complete derivation of the eigenoperators of the linearized 
approximation to the renormalization group transformation about the critical Gaussian 
fixed point, but the derivation is simple. 4 Increasing powers of transverse derivatives and 
increasing powers of creation and annihilation operators lead to increasingly irrelevant 
operators. The irrelevant operators are called 'non-renormalizable' in old-fashioned 
Feynman perturbation theory. Their magnitude decreases at an exponential rate as 
the cutoff is lowered, which means that they increase at an exponential rate as the 
cutoff is raised and produce increasingly large divergences if we try to follow their 
evolution perturbatively in this exponentially unstable direction. 

The only relevant operator is the mass operator, 

fi 2 J dq, dq 2 (167r 3 )<5 3 ( gi - q 2 ) a\ qi )a(q 2 ) , (69) 

while the fixed point Hamiltonian is marginal (of course), and the operator in which 
U4 = A (a constant) is marginal. A 3 operator would also be relevant. 

The next logical step in a renormalization group analysis is to study the transfor- 
mation to second order in the interaction, keeping the second-order corrections from 
the similarity transformation. I will compute the correction to u 2 to this order and 
refer the interested reader to Ref. 4 for more complicated examples. 

The matrix element of the one-body operator between single particle states is 

(p'\h\p) = (0| a{p') h a\p) |0> = (167T 3 )5V - p) u 2 (-p,p) ■ (70) 

Thus, we easily determine u 2 from the matrix element. It is easy to compute matrix 
elements between other states. We computed the matrix elements of the effective 
Hamiltonian generated by the similarity transformation when the cutoff is lowered in 
Eq. (57), and now we want to compute the second-order term generated by the four- 
point interactions above. There are additional corrections to u 2 at second order in the 
interaction if Uq, etc. are nonzero. 

Before rescaling we find that the transformed Hamiltonian contains 

(lQir 3 )8 3 (p' -p) 5v 2 (-p,p) = 



hS dk 



dfadfa 9(A — \p — k 1 — k 2 — k 3 |) 



e(lP - _ k: _ k - _ ,-| _ , Ao) VWK *..*.><*■■ *>Mm , (71) 

p k 1 k 2 k 3 
where p~ = Pj_/p + , etc. One can readily verify that 

(p'\V\fa, k 2 , k 3 ) = (lQn 3 )5 3 (p' -h-k 2 - k 3 ) e(p~ - fcr - k 2 - fcg ) u 4 (-p', h, k 2 , k 3 ) , 

(72) 

(fa, fa, h\V\p) = (16n 3 )5 3 (p-fa-k 2 -k 3 ) Q(p -fa -k 2 -k 3 ) u 4 (-fa, -k 2 , -k 3 ,p) . 

(73) 

This leads to the result, 

8v 2 (-p,p) = J dfadk 2 dk 3 (\Qti 3 )5 3 (p — fa—k 2 — k 3 ) 

6(A - \p~ - fcf - lq - fcg I) G(\p- - fcf -lq r]A ) 

M-p, fa, fa, fa) M-fa, -fa, -fa,p) ^ 74 ^ 

p — fa — fa — k 3 
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To obtain 5u2 from 5v2 we must rescale the momenta, the fields, and the Hamiltonian. 
The final result is 




(75) 



where 



P ~ k-i — k-2 — k% 
p+' = p+ ; p 1 -' = rjp- 1 , kf = kf , and kf = rjkf . 



Coupling Coherence 

The basic mathematical idea behind coupling coherence was first formulated by 
Oehme, Sibold, and Zimmerman. 33 They were interested in field theories where many 
couplings appear, such as the standard model, and wanted to find some means of 
reducing the number of couplings. Wilson and I developed the ideas independently in 
an attempt to deal with the functions that appear in marginal and relevant light-front 
operators. 3 

The puzzle is how to reconcile our knowledge from covariant formulations of QCD 
that only one running coupling constant characterizes the renormalized theory with 
the appearance of new counterterms and functions required by the light-front formu- 
lation. What happens in perturbation theory when there are effectively an infinite 
number of relevant and marginal operators? In particular, does the solution of the per- 
turbative renormalization group equations require an infinite number of independent 
counterterms {i.e., independent functions of the cutoff)? Coupling coherence provides 
the conditions under which a finite number of running variables determines the renor- 
malization group trajectory of the renormalized Hamiltonian. To leading nontrivial 
orders these conditions are satisfied by the counterterms introduced to restore Lorentz 
covariance in scalar field theory and gauge invariance in light-front gauge theories. In 
fact, the conditions can be used to determine all counterterms in the Hamiltonian, 
including relevant and marginal operators that contain functions of longitudinal mo- 
mentum fractions; and with no direct reference to Lorentz covariance, this symmetry 
is restored to observables by the resultant counterterms in scalar field theory. 4 

A coupling-coherent Hamiltonian is analogous to a fixed point Hamiltonian, but 
instead of reproducing itself exactly it reproduces itself in form with a limited number 
of independent running couplings. If qa is the only independent coupling in a theory, 
in a coupling-coherent Hamiltonian all other couplings are invariant functions of g^, 
fi(g\). The couplings fi(g\) depend on the cutoff only through their dependence on 
the running coupling g^, and in general we demand /i(0) = 0. This boundary condition 
on the dependent couplings is motivated in our calculations by the fact that it is the 
combination of the cutoff and the interactions that force us to add the counterterms 
we seek, so the counterterms should vanish when the interactions are turned off. 

Let me start with a simple example in which there is a finite number of relevant 
and marginal operators ab initio, and use coupling coherence to discover when only one 
or two of these may independently run with the cutoff. In general such conditions are 
met only when an underlying symmetry exists. 

Consider a theory in which two scalar fields interact, 



(76) 
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Under what conditions will there be fewer than three independent running coupling 
constants? We can use a simple cutoff on Euclidean momenta, q 2 < A 2 . Letting 
t = ln(A/A ), the Gell-Mann-Low equations are 

f)\ 1 

-Of = 3CA? + -CA^ + 0(2 loop) , (77) 
f)\ 1 

-gf = 3CA^ + -CA^ + 0(2 loop) , (78) 
<9A 2 

= gCA§ + CAiA 3 + CA 2 A 3 + 0(2 loop) ; (79) 

where £ = h/ (16-7T 2 ). It is not important at this point to understand how these equations 
are derived. 

First suppose that Ai and A2 run separately, and ask whether it is possible to find 
As(Ai, A 2 ) that solves Eq. (79). To one-loop order this leads to 

( 3A ? + ^) It + ( 3A * + is A ») k = + AiA > + A ^ • (80) 

If Ai and A2 are independent, we can equate powers of these variables on each side of 
Eq. (80). If we allow the expansion of A3 to begin with a constant, we find a solution 
to Eq. (80) in which all powers of Ai and A2 appear. In this case a constant appears on 
the right-hand-sides of Eqs. (77) and (78), and there will be no Gaussian fixed points 
for Ai and A 2 . We are generally not interested in the possibility that a counterterm 
does not vanish when the canonical coupling vanishes, so we will simply discard this 
solution both here and below. We are interested in the conditions under which one 
variable ceases to be independent, and the appearance of such an arbitrary constant 
indicates that the variable remains independent even though its dependence on the 
cutoff is being reparameterized in terms of other variables. 

If we do not allow a constant in the solution, we find that A3 = a\\ + /3A 2 + 0(X 2 ). 
When we insert this in Eq. (80) and equate powers on each side, we obtain three 
coupled equations for a and j3. These equations have no solution other than a = and 
f3 = 0, so we conclude that if Ai and A 2 are independent functions of t, A 3 will also be 
an independent function of t unless the two fields decouple. 

Assume that there is only one independent variable, A = Ai, so that A 2 and A3 are 
functions of A. In this case we obtain two coupled equations, 

( 3 A 2 + 1a 2 ) ®h = 3A 2 + —A 2 , (81) 

(3A 2 + ^A 2 ) d -h = ^A 2 + AA 3 + A 2 A 3 . (82) 

If we again exclude a constant term in the expansions of A 2 and A3 we find that the 
only non-trivial solutions to leading order are A 2 = A, and either A3 = 2A or A3 = 6A. 
If A 3 = 2A, 

V(4>) = ^(<Pl + rtf , (83) 
and we find the 0(2) symmetric theory. If A3 = 6 A, 



2^4! 



v (<t>) = 7T-n + + (0i - <^) 4 1 , (84) 
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and we find two decoupled scalar fields. Therefore, A2 and A3 do not run independently 
with the cutoff if there is a symmetry that relates their strength to Ai. 

The condition that a limited number of variables run with the cutoff does not only 
reveal symmetries broken by the regulator, it may also be used to uncover symmetries 
that are broken by the vacuum. I will not go through the details, but it is straight- 
forward to show that in a scalar theory with a 3 coupling, this coupling can be fixed 
as a function of the <fi 2 and 4 couplings only if the symmetry is spontaneously broken 
rather than explicitly broken. 3 

This example is of some interest in light-front field theory, because it is difficult 
to reconcile vacuum symmetry breaking with the requirements that we work with a 
trivial vacuum and drop zero-modes in practical non-perturbative Hamiltonian calcu- 
lations. Of course, the only way that we can build vacuum symmetry breaking into the 
theory without including a nontrivial vacuum as part of the state vectors is to include 
symmetry breaking interactions in the Hamiltonian and work in the hidden symme- 
try phase. The problem then becomes one of finding all necessary operators without 
sacrificing predictive power. The renormalization group specifies what operators are 
relevant, marginal, and irrelevant; and coupling coherence provides one way to fix the 
strength of the symmetry-breaking interactions in terms of the symmetry-preserving 
interactions. This does not solve the problem of how to treat the vacuum in light-front 
QCD by any means, because we have only studied perturbation theory; but this result 
is encouraging and should motivate further investigation. 

For the QED and QCD calculations discussed below, I need to compute the hamil- 
tonian to second order, while the canonical coupling runs at third order. To determine 
the generic solution to this problem, I present an oversimplified analysis in which there 
are three coupled renormalization group equations for the independent marginal cou- 
pling (g), in addition to dependent relevant (/i) and irrelevant (w) couplings. 

A*m«+i = 4 /i m! + Cn9i + 0(gf) , (85) 

9m l+1 = 9n H + c g gf + 0{gf) , (86) 

w mi+1 = \w mi + c w gf + O(gf) . (87) 

I assume that /ii = (gf + O(gf) and wi = rjgf + O(gf), which satisfy the conditions 
of coupling coherence. Substituting into the renormalization group equations yields 
(dropping all terms of O(gf)), 

Cgf+i = (a? = Kgl + , (88) 
wf+i = vgf = \vgf + c w gf . (89) 

The solutions are ( = — |c M and rj = |c w . 

These are exactly the coefficients in a Taylor series expansion for /i(g) and w(g) 
that reproduce themselves. This observation suggests an alternative way to find the 
coupling coherent Hamiltonian without explicitly setting up the renormalization group 
equations. Although this method is less general, we only need to find what operators 
must be added to the Hamiltonian so that at 0(g 2 ) it reproduces itself, with the only 
change being the change in the specific value of the cutoff. Coupling coherence allows 
us to substitute the running coupling in this solution, but it is not until third order 
that we would explicitly see the coupling run. This is how I will fix the QED and QCD 
hamiltonians to second order. 
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QED and QCD Hamiltonians 

In order to derive the renormalized QED and QCD Hamiltonians I must first list 
the canonical Hamiltonians. I follow the conventions of Brodsky and Lepage. 34 There 
is no need to be overly rigorous, because coupling coherence will fix any perturbative 
errors. I recommend the papers of Zhang and Harindranath 6 for a more detailed dis- 
cussion from a different point of view. The reader who is not yet concerned with details 
can skip this section. 

I will use A + = gauge, and I drop zero modes. I use the Bjorken and Drell 
conventions for gamma matrices. The gamma matrices are 

where are the Pauli matrices. This leads to 

H-la Hi -?)■ < 9i > 

Useful identities for many calculations are 7+7~7 + = 47 + , and 7^7+7^ = 47^. 
The operator that projects onto the dynamical fermion degree of freedom is 



1 + 1 - + 1 ( 1 
A +=2 7V = 4 7 7 = 2 U 



03 
<7 3 1 



(92) 



and the complement projection operator is 



and, 



The Dirac spinors u{p, a) and v(p, a) satisfy 

(/ - m)u(p, a) — , + m)v(p, a) — , (94) 

o-)m(p, ct') = -v(p, a)v(p, a') = 2m5 aa > , (95) 
u(p, (T)fu(p, a') = u(p, a^vip, O = 2^5^/ , (96) 
u(p,a)u(p,a) =tf + m , ^ v(p,<j)v(p,a) — m . (97) 

There are only two physical gluon (photon) polarization vectors, eu_ and e 2 ±] but 
it is sometimes convenient (and dangerous once covariance and gauge invariance are 
violated) to use e M , where 

e + = 0, e~ = ^^. (98) 

It is often possible to avoid using an explicit representation for e±, but completeness 
relations are required, 

E£(A)tf(A) = -sr, (99) 

A 

so that, 

Ee"(A)e*"(A) = -«?f + + «) + » ( 10 °) 



A 
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where i] + = rj 1 = rf = and i]~ = 2. One often encounters diagrammatic rules in 
which the gauge propagator is written so that it looks covariant; but this is dangerous 
in loop calculations because such expressions require one to add and subtract terms 
that contain severe infrared divergences. 
The QCD Lagrangian density is 



(101) 



where F» v = d»A v - d v A> x + ig[A», A u ] and iD* = id» - gA^. The SU(3) gauge fields 
are A^ = I] a A^T a , where T a are one-half the Gell-Mann matrices, A a , and satisfy 
Tr T a T b = 1/2 5 ab and [T a ,T b ] = if abc T c . 

The dynamical fermion degree of freedom is ip + = A + ip, and this can be expanded 
in terms of plane wave creation and annihilation operators at x + = 0, 



dk+d 2 kj 



k+>o 167r 3 A; + 



b r (k, a)u+(k, a)e~ ik ' x + d r \k, a)v + (k, a)e 



„ik-x 



<r=±l/2 ' 

where these field operators satisfy 

{V + {x)^{y)} x+=y+=Q = A + S rs 6 3 (x - y) , 
and the creation and annihilation operators satisfy 

{b r (k, a), 6 at (p,</)} = {d r (k,a),d^(p,a')} = 16n 3 k + 5 r J aa ,6 3 (k - p) . 



(102) 



(103) 



(104) 



The indices r and s refer to SU(3) color. In general, when momenta are listed without 
specification of components, as in S 3 (p), I am referring to p + and p_i_. 

The transverse dynamical gluon field components can also be expanded in terms 
of plane wave creation and annihilation operators, 



A 



dk + d 2 k± 

k+>0 167T 3 fc + 



a c (k, \)e\(\)e- ik - x + a ct (/c, A)e i _ L *(A)e 



Ak-x 



(105) 



The superscript % refers to the transverse dimensions x and y, and the superscript c is 
for SU(3) color. If required the physical polarization vector can be represented 



cl(t) = -^(M) , c±a) = ^(i,-o- 



The quantization conditions are 

'A i l{x),d + A j l{y) 



= i8 ij S cd S 3 (x - y) , 



x+=y+=0 

a c (k,X),a d \p,X')} = 167i 3 k + 5 X y5 cd 5 3 (k-p) 



(106) 



(107) 



(108) 



The classical equations for vf)_ = A_ip and A~ do not involve time-derivatives, 
so these variables can be eliminated in favor of dynamical degrees of freedom. This 
formally yields 



1 



id+ 



ia± ■ D_|_ + [3m ip^ 



(109) 
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where the variable ■?/>_ is defined on the second line to separate the interaction-depen- 
dent part of ■?/>_; and 

= ~A a ~ + M^^id + A b l)A c l + 2^T>+ j , (110) 

where the variable A~ is defined on the second line to separate the interaction-depen- 
dent part of A~ . 

Given these replacements, we can follow a canonical procedure to determine the 
Hamiltonian. This path is full of difficulties that I ignore, because ultimately I will 
use coupling coherence to refine the definition of the Hamiltonian and determine the 
non-canonical interactions that are inevitably produced by the violation of explicit 
covariance and gauge invariance. For my purposes it is sufficient to write down a 
Hamiltonian that can serve as a starting point: 

H = H + V, (111) 



H = J ^x^ridiAidiA^) + ^{(ia l ± d l ± + /3m)^+} 
dk + d 2 k ± 



= E 

colors 



167r 3 Ar 



£^at(MMM) 



+ E k ^ m2 ( 6t ( fc » a ) b (k, <?) + d\k, a)d(k, a 



;ii2) 



In the last line the 'self-induced inertias' (i.e., one-body operators produced by normal- 
ordering V) are not included. It is difficult to regulate the field contraction encountered 
when normal-ordering in a manner exactly consistent with the cutoff regulation of 
contractions encountered later. Coupling coherence avoids this issue and produces the 
correct one-body counterterms with no discussion of normal-ordering required. 

The interactions are complicated and are most easily written using the variables, 
■0 = tjj- + and A, where A + = 0, A~ is defined above, and A l ± = A\. Using these 
variables we have 



V = d A x 



gfapA^ + 2g Tr(id^[A^A u ]) 



— 1 r 
2 



A»A V 



1 



A^A V )+ g ^A^^_^A^ 



2^,+ i 



-g ^7 



1 



(id+y 



+g 2 Tri id + A^A^ 



(id+Y 



id + A u , A v 



(113) 



The commutators in this expression are SU(3) commutators only. The potential al- 
gebraic complexity of calculations becomes apparent when one systematically expands 
every term in V and replaces: 



4>- 



id+ 



(114) 
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A-^T^id ± -A ± ; (115) 

and then expands ip + and Aj_ in terms of creation and annihilation operators. It rapidly 
becomes evident that one should avoid such explicit expansions if possible. 



LIGHT-FRONT QED 

In this section I will follow the strategy outlined in the first section to compute 
the positronium spectrum. I will detail the calculation through the leading order Bohr 
results 8 and indicate how higher order calculations proceed. 10 ' 11 

The first step is to compute a renormalized cutoff Hamiltonian as a power series in 
the coupling e. Starting with the canonical Hamiltonian as a 'seed,' this is done with 
the similarity renormalization group 17, 18 and coupling coherence. 3 ' 4 The result is an 
apparently unique perturbative series, 

= h + e A hx + e\h 2 + • • • + efh N . (116) 

Here e A is the running coupling constant, and all remaining dependence on A in the 
operators hi must be explicit. In principle I must also treat m\, the electron running 
mass, as an independent function of A; but this will not affect the results to the order I 
compute here. We must calculate the Hamiltonian to a fixed order, and systematically 
improve the calculation later by including higher order terms. 

Having obtained the Hamiltonian to some order in e, the next step is to split it 
into two parts, 

H A = H A + V A . (117) 

As discussed before, Ti A must be accurately solved non-perturbatively, producing a ze- 
roth order approximation for the eigenvalues and eigenstates. The greatest ambiguities 
in the calculation appear in the choice of , which requires one of science's most 
powerful computational tools, trial and error. 

In QED and QCD I assume that all interactions in H A preserve particle number, 
with all interactions that involve particle creation and annihilation in V A . This assump- 
tion is consistent with the original hypothesis that a constituent picture will emerge, 
but it should emerge as a valid approximation. 

The final step before the loop is repeated, starting with a more accurate approxi- 
mation for H A , is to compute corrections from V A in bound state perturbation theory. 
There is no reason to compute these corrections to arbitrarily high order, because the 
initial Hamiltonian contains errors that limit the accuracy we can obtain in bound state 
perturbation theory. 

In this section I: (i) compute H A to 0(e 2 ), (ii) assume the cutoff is in the range 
am 2 < A 2 < a 2 m 2 for non-perturbative analyses, (iii) include the most infrared singular 
two-body interactions in TLq, and (iv) estimate the binding energy for positronium to 
0{a 2 m). 

Since H A is assumed to include interactions that preserve particle number, the 
zeroth order positronium ground state will be a pure electron-positron state. We only 
need one- and two-body interactions; i.e., the electron self-energy and the electron- 
positron interaction. The canonical interactions can be found in Eq. (113), and the 
second-order change in the Hamiltonian is given in Eq. (57). The shift due to the bare 
electron mixing with electron-photon states to lowest order (see Figure 1) is 

2 , dkfd 2 k 1± 6{k+)6{kt) u^D^k^rifc + m)Yu(p) 

02jp — 6 



J 16tt 3 kfkt p'-ki-fc 
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where, 



p + - . 

for = 



**2_L 



p ± - k 



1± 



fof ± + m? 



(k+) 



A, 



+ 



+ p: (fo±^ + fc^) , 



(119) 
(120) 

(121) 
(122) 



I have not yet displayed the cutoffs. To evaluate the integrals it is easiest to use Jacobi 
variables x and s for the relative electron-photon motion, 



fo+ = X p + , ki_L = xp± + s , 



which implies 



^2 



[l-x)p + , k 2 _L = (1 - x)p ± - s . 
The second-order change in the electron self-energy becomes 



(123) 
(124) 



dxd 2 s 
16tt 3 
1 



J a 2 s 2 + x 2 m 2 \ n ( s 2 + x 2 m 2 
9[yA 2 -|0' 



x(l — x) 



s 2 + x 2 m 2 



x(l — x) 
u(p, cr) <^2(1 — x) p 1 — 2m + 

./+ 

x 



-yA? 



2_ 



2^2 

,2 



2s 2 x(2 — x)m 2 



>u(p,a) 



(125) 



where y = p + / P + . 

It is straightforward in this case to determine the self-energy required by coupling 
coherence. Since the electron-photon coupling does not run until third order, to second 
order the self-energy must exactly reproduce itself with A — > For the self-energy 
to be finite we must assume that ST, reduces a positive self-energy, so that 



e£*(p) 



e\ r dxd 2 s 



16tt 3 



9(yx-e)9[yA 2 - 



u(p,a)l2(l-x) i>-2m + ^- 
I ' p+ 



s 2 + x 2 m 2 \ 
x(l — x) J 
2s 2 Js 2 

2~ " 



s 2 + x 2 m 2 J 
x(2 — x)m 2 



x. 



y 2 A 2 



8?rV 



2yA 2 In! 



m?)e 



(yA 2 + 
+ 3m 2 In 



yA 



yA 2 + m 2 

I have been forced to introduce a second cutoff, 

xp + > eP + , 



1 — x 1 — X 

2 1 ym 2 A 2 
2 yA 2 + m 2 

'\+0(e/y). 



>u(p, a) 



(126) 



(127) 



because after the s integration is completed we are left with a logarithmically divergent 
x integration. Other choices for this second infrared cutoff are possible and lead to 
similar results. This second cutoff must be taken to zero and no new counterterms can 
be added to the Hamiltonian, so all divergences must cancel before it is taken to zero. 

The electron and photon (quark and gluon) 'mass' operators, are a function of 
a longitudinal momentum scale introduced by the cutoff, and there is an exact scale 
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invariance required by longitudinal boost invariance. Here I mean by 'mass operator' 
the one-body operator when the transverse momentum is zero, even though this does 
not agree with the free mass operator because it includes longitudinal momentum de- 
pendence. The cutoff violates boost invariance and the mass operator is required to 
restore this symmetry. 

We must interpret this new infrared divergence, because we have no choice about 
whether it is in the Hamiltonian if we use coupling coherence. We can only choose 
between putting the divergent operator in TCq or in V A . I make different choices in 
QED and QCD, and the arguments are based on physics. 

The divergent electron 'mass' is a complete lie. We encounter a term proportional 
to e\A 2 ln(l/e)/P + when the scale is A; however, we can reduce this scale as far as 
we please in perturbation theory. Photons are massless, so the electron will continue 
to dress itself with small-x photons to arbitrarily small A. Since I believe that this 
divergent self-energy is exactly canceled by mixing with small-x photons, and that 
this mixing can be treated perturbatively in QED, I simply put the divergent electron 
self-energy in V A , which is treated perturbatively. 

There are two time-ordered diagrams involving photon exchange between an elec- 
tron with initial momentum p± and final momentum p 2 , and a positron with initial 
momentum k\ and final momentum k 2 . These are shown in Figure 4, along with the 
instantaneous exchange diagram. Using Eq. (57), we find the required matrix element 
of SH, 

e 2 

SH = -—D flI/ (q)u(p 2 , a 2 )Yu(Pi, cr i)v(kiAih u v(k 2 , X 2 ) 




0(\pI -P2-q~\- a 2 /^ + ) e(\ P j - p 2 - q -\ - \k 2 - fcr - q -\) 

Pi - P2 - <r 

0(\k 2 ~K-q-\- Aj/V + ) 0(|A; 2 - -K-q~\- br - P2 - q~\) 



+ 



k 2 — ki — q 

o(§^-\Pi+K-P2-k 2 |) , (128) 



where q + = pf - p 2 , q± = p 1± - p 2 ±, and q = q\/q + . 

1 have used the second cutoff on longitudinal momentum that 1 was forced to 
introduce when computing the change in the self-energy. We will see in the section on 
confinement that it is essential to include this cutoff everywhere consistently. In QED 
this point is not immediately important, because all infrared singular interactions, 
including the infrared divergent self-energy, are put in V A and treated perturbatively. 
Divergences from higher orders in V A cancel. 10, 11 

To determine the interaction that must be added to the Hamiltonian to maintain 
coupling coherence, we must again find an interaction that when added to SV repro- 
duces itself with Ao — > Ai everywhere. The coupling coherent interaction generated by 
the first terms in SH are not uniquely determined at this order. There is some ambigu- 
ity because we can obtain coupling coherence either by having SH increase the strength 
of an operator by adding additional phase space strength, or we can have SH reduce the 
strength of an operator by subtracting phase space strength. The ambiguity is resolved 
in higher orders, so I will simply state the result. If an instantaneous photon-exchange 
interaction is present in H, SH cancels part of this marginal operator and increases the 
strength of a new photon-exchange interaction. This new interaction reproduces the 
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i q 



Figure 4. Effective two-body interaction between low energy constituents resulting from: (i) the 
canonical instantaneous exchange interaction, and (ii) the elimination of direct coupling between low 
energy two-body states and high energy states containing an additional gauge particle. 



effects of high energy photon exchange removed by the cutoff. The result is 



^coh 



^D llv {q)u(p 2 , a 2 )Yu(pi, cr^h, \ l )Yv{k 2 , A 2 ) 

e{\ q + \ - ep+) e(^- 1 P r + K -P2 - K |) 

0(\Pi -P2-q~\- A 2 /^ + ) 0{\Pi -P2-q~\- \K ~ K - q~\) 



Pi -P2 - or 



+ 



QJ\K - K - gl - A 2 /^ + ) 9(\K -K-q-\- IpT - vi - fl) ' 

k 2 — ki — q 



(129) 



This matrix element exactly reproduces photon exchange above the cutoff. The 
cutoff removes the direct coupling of electron-positron states to electron-positron-pho- 
ton states whose energy differs by more than the cutoff, and coupling coherence dictates 
that the result of this mixing should be replaced by a direct interaction between the 
electron and positron. We could obtain this result by much simpler means at this order 
by simply demanding that the Hamiltonian produce the 'correct' scattering amplitude 
at C(e 2 ) with the cutoffs in place. Of course, this procedure requires us to provide the 
'correct' amplitude, but this is easily done in perturbation theory. 

V^ h is non-canonical, and we will see that it is responsible for producing the 
Coulomb interaction. We need some guidance to decide which irrelevant operators are 
most important. We find a posteriori that differences of external transverse momenta, 
and differences of external longitudinal momenta are both proportional to a. This 
allows us to identify the dominant operators by expanding in powers of these implicit 
powers of a. This indicates that it is the most infrared singular part of V co h that is 
important. As explained above, this operator receives substantial strength only from 
the exchange of photons with small longitudinal momentum; so we expect inverse q + 
dependence to indicate 'strong' interactions between low energy pairs. So the part of 
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V co h that is included in Hq is 



v 



q 2 

e ^ (^+)3 a ih +u (P^ ff iM fc i> ^ih + v(k 2 , A 2 ) 



co/i 



/A 2 

eP+) ^IpT" |Pl +kl ~ P2 ~ k * 



0(\PI -P2-Q-\- A 2 /^ + ) ^(br -P2-Q-\- \K ~ K ~ q~\) 



Pi -P2 -<r 



+ 



e(\k£ - k - <r\ - a 2 /p + ) g(|fcj -K- g -\- br - P2 - q 

^2 ~ ^1 ~~ Q 



/A 2 



it? 

e(\q + \-eP + ) e 



Pi +hi -p 2 -k 2 



0{\Pi -P2-q~\- A 2 /V + ) 9(\pT -P2-q~\- \k 2 - K ~ 



+ 



\K -k - q -\-A 2 /r- 



Pi -P2 -or 



h ~q | - \Pi ~P 2 ~q I) 



k-2 ~ k\ — q 
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The Hamiltonian is almost complete to second order in the electron-positron sector, 
and only the instantaneous photon exchange interaction must be added. The matrix 
element of this interaction is 



Vi ns tant = SU — u(p 2 , <7 2 )7 + u(pi , G\ )v{h , Ai )l + v(k 2 , A 2 ) 



/ A 2 

x 6 [^T- \Pi +k x -p 2 -k 2 



= -^A^PlPtkt kt S <Tl<T2 S \l 



A 2 



p 1 +k 1 -p 2 - k 2 



(131) 



The only cutoff that appears is the cutoff directly run by the similarity transformation 
that prevents the initial and final states from differing in energy by more than A 2 /P + . 

This brings us to a final subtle point. Since there are no cutoffs in Vi nsta nt that di- 
rectly limit the momentum exchange, the matrix element diverges as q + — > 0. Consider 
V co h in this same limit, 



V, 



coh 



4e A Jpip 2 ktkJ — 5 aia2 S Xl x 2 



X 9 | q 



A 2 ^ 
P+ 



'A*_ 

p+ 



Pi + h - p 2 - k 2 



(132) 



This means that as q + — > 0, V co h partially screens V instant , leaving the original operator 
multiplied by 9(A 2 /P + — \q~\). However, even after this partial screening, the matrix 
elements of the remaining part of Vinstant between bound states diverge and we must 
introduce the same infrared cutoff used for the self-energy to regulate these divergences. 
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This is explicitly shown in the section on confinement. However, all divergences from 
Vinstant are exactly canceled by the exchange of massless photons, which persists to 
arbitrarily small cutoff. This cancellation is exactly analogous to the cancellation of 
the infrared divergence of the self-energy, and will be treated in the same way. The 
portion of V instant that is not canceled by V co h will be included in V A , the perturbative 
part of the Hamiltonian. We will not encounter this interaction until we also include 
photon exchange below the cutoff perturbatively, so all infrared divergences should 
cancel in this bound state perturbation theory. I repeat that this is not guaranteed for 
arbitrary choices of 7i A , and we are not free to simply cancel these divergent interactions 
with counterterms because coupling coherence completely determines the Hamiltonian. 

We now have the complete interaction that I include in 7Y A . Letting 7i A = h + Vq, 
where ho is the free hamiltonian, I add parts of Vi nsta nt and V co h to obtain 

/ / A 2 




v o = - 4e A\/PiP2^i kj 5 aia2 5 Xl x 2 9\ — - | Pi +k ± -p 2 - k 2 | ] 



+ 



0(\Pi-p 2 -q-\-A 2 /V 



+ 



(q + ) 3 Pi -P2 -q (q + ) 2 _ 

x 0(\pT -P2 ~Q~\ - \ k 2 ~K -Q~\) 
n ' + T ^\e(\k 2 -K-q-\-K 2 /V + ) 



(q + ) 3 K-K-q- (q+)\ 

X 6(\k 2 - fcf - q-\ - fa - p 2 - . (133) 

In order to present an analytic analysis I will make assumptions that can be jus- 
tified a posteriori. First I will assume that the electron and positron momenta can be 
arbitrarily large, but that in low-lying states their relative momenta satisfy 

|p± - k ± | < m , (134) 

\p + -k + \ <^p + + k + . (135) 

It is essential that the condition for longitudinal momenta not involve the electron mass, 
because masses have the scaling dimensions of transverse momenta and not longitudinal 
momenta. As above, I use p for the electron momenta and k for the positron momenta. 
To be even more specific, I will assume that 

|p_i_ — kj_| ~ am , (136) 

\p + -k + \ ~a(p + + k + ) . (137) 

This allows us to use power counting to evaluate the perturbative strength of operators 
for small coupling, which may prove essential in the analysis of QCD. 1 Note that these 
conditions allow us to infer 

|Pi_l - P2_l| ~ am , (138) 

\pf~Pt\ ~a(p + + k + ) . (139) 

Given these order of magnitude estimates for momenta, we can drastically simplify 
the free energies in the kinetic energy operator and the energy denominators in Vq . We 
can use transverse boost invariance to choose a frame in which 

Pi = Vi p+ , Pi± = «i , K = (1 - Vi)P + , k i± = -Ki , (140) 
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so that 



P + (Pi -P 2 ~Q ) = 



k\ + m 2 



k\ + m 2 



[«1 - K 2 



yi 

= -Am 2 {yi 



yi 



2/2) 



2/i — 3/2 

2 



+ 0(o; 2 m 2 ) . 



(141) 



2/1 - y2 

To leading order all energy denominators are the same. Each energy denominator 
is 0(am 2 ), which is large in comparison to the binding energy we will find. This is 
important, because the bulk of the photon exchange that is important for the low energy 
bound state involves intermediate states that have larger energy than the differences in 
constituent energies in the region of phase space where the wave function receives most 
of its strength. This allows us to use a perturbative renormalization group to compute 
the dominant effective interactions. 

There are similar simplifications for all energy denominators. After making these 
approximations we find that the matrix element of Vq is 

\2 

A 2 ) 



m - 2/2I + 



[«1 - K 2 



v v lS/i-2/2| 

e^K 2 - 4 I k\ + 4m 2 (l - 2 yi f -k\- 4m 2 (l - 2y 2 f 

(«i - r 1 



= -\Qe\m 2 \/ yiy 2 {l 



(yi - y 2 ) 2 [4m 2 (?/i - y 2 ) 2 + («i - k 2 ) 2 (/ci - k 2 ) 2 

[ K l ~ A 2 



2/2 1 



yi)(l-y 2 ) ^(4m 2 || /l 

fl(A 2 - 4 I k\ + 4m 2 (l - 2 yi f -k\- 4m 2 (l - 2y 2 ) 

1 



Vl-V2\ 

2 
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_4m 2 (y 1 - y 2 ) 2 + («i - k 2 ) 2 _ 

In principle the electron-positron annihilation graphs should also be included at this 
order, but the resultant effective interactions do not diverge as q + — > 0, so I include 
such effects perturbatively in V A . 

At this point we can complete the zeroth order analysis of positronium using the 

state, 

<f>(p,a;k,Xp(p,a)S(k,X)\0) , (143) 

where <f>(p, a; k, A) is the wave function for the relative motion of the electron and 
positron, with the center-of-mass momentum being P. We need to choose the longi- 
tudinal momentum appearing in the cutoff, and I will use the natural scale P + . The 
matrix element of TCq is 

(tf (P)|7#|tt(P')) = 16tt 3 5 3 (P - P') x 
dyd 2 K 



16tt 3 



4m 2 + 4k 2 + 4m 2 (l - 2y) 2 \<f>{K,y)\ 



16e 2 m 2 



dy\d 2 K\ dy 2 d 2 K 2 

167T 3 167T 3 



9 Um 2 \ yi -y 2 \ + 



[Kl ~ K 2 ) 

\yi -2/2 1 

#(A 2 -A\k\ + 4m 2 (l - 2 yi f 4m 2 (l - 2y 2 )' 

1 



A 2 ) 



4m 2 |y 1 - y 2 \ 2 + («i - k 2 ) 2 



0*(k2,S/2)0(ki,2/i) 



(144) 
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I have chosen a frame in which P± = and used the Jacobi coordinates defined above, 
and indicated only the electron momentum in the wave function since momentum 
conservation fixes the positron momentum. I have also dropped the spin indices because 
the interaction in Hq is independent of spin. 

If we vary this expectation value subject to the constraint that the wave function 
is normalized we obtain the equation of motion, 

M 2 0( Kl , yi) = {Am 2 - AmE + £ 2 )0(«i, Vl ) 

= Am 2 + AK 2 1 +Am 2 (l-2y 1 ) 2 ](p( y K 1 ,y 1 ) 

-16e m / 9 [Am \y x - y 2 \ + -, r - A 

Q{k 2 - A | k\ + 4m 2 (l - 2 Vl ) 2 - k\ - Am 2 (\ - 2y 2 ) 2 |) 

-j— 2i ^—7 ^0(« 2 ,y 2 )- (145) 

_Am 2 \y 1 - y 2 \ 2 + («i - /t 2 ) 2 J 

S is the binding energy, and we can drop the -E 2 term since it will be 0(a 4 ). 

I do not think that it is possible to solve this equation analytically with the cutoffs 
in place, and with the light-front kinematic constraints — 1 < 1 — 2y, < 1. In order to 
determine the binding energy to leading order, we need to evaluate the regions of phase 
space removed by the cutoffs. 

If we want to find a cutoff for which the ground state is dominated by the electron- 
positron component of the wave function, we need the first cutoff to remove the 
important part of the electron-positron-photon phase space. Using the 'guess' that 
\k\ = 0(am) and 1 — 2y = 0(a), this requires 

A 2 < am 2 . (146) 

On the other hand, we cannot allow the cutoff to remove the region of the electron- 
positron phase space from which the wave function receives most of its strength. This 
requires 

A 2 > a 2 m 2 . (147) 

While it is not necessary, the most elegant way to proceed is to introduce 'new' 
variables, 

Ki = k ±i , (148) 

y t = -+ - h . (149) 

2 2,/k 2 + jfc2 + m 2 



This change of variables can be 'discovered' in a number of ways, but they basically 
take us back to equal time coordinates, in which both boost and rotational symmetries 
are kinematic after a nonrelativistic reduction. 

For cutoffs that satisfy am 2 > A 2 > a 2 m 2 , Eq. (145) simplifies tremendously 
when all terms of higher order than a 2 are dropped. Using the scaling behavior of the 
momenta, and the fact that we will find E 2 is 0(a 4 ), Eq. (145) reduces to: 

k 2 r rl 3 k 1 

-*i^«-»/(4(M* w ' (150) 

The step function cutoffs drop out to leading order, leaving us with the familiar non- 
relativistic Schrodinger equation for positronium in momentum space. The solution 
is 

m = 77 2 f F , 2 , (i5i) 

(k 2 + mE) 
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E = -a 2 m . (152) 

M is a normalization constant. 

This is the Bohr energy for the ground state of positronium, and it is obvious that 
the entire nonrelativistic spectrum is reproduced to leading order. 

Beyond this leading order result the calculations become much more interesting, 
and in any Hamiltonian formulation they rapidly become complicated. There is no 
analytic expansion of the binding energy in powers of a, since negative mass-squared 
states appear to signal vacuum instability when a is negative ; but our simple strategy 
for performing field theory calculations can be improved by taking advantage of the 
weak coupling expansion. 1 

We can expand the binding energy in powers of a, and as is well known we find 
that powers of log (a) appear in the expansion at 0(a 5 ). We have taken the first step to 
generate this expansion by expanding the effective Hamiltonian in the explicit powers 
of a which appear in the renormalization group analysis. The next step is to take ad- 
vantage of the fact that all bound state momenta are proportional to a, which allows us 
to expand each of the operators in the Hamiltonian in powers of momenta. The renor- 
malization group analysis justifies an expansion in powers of transverse momenta, and 
the nonrelativistic reduction leads to an expansion in powers of longitudinal momentum 
differences. The final step is to regroup terms appearing in bound state perturbation 
theory. For example, when we compute the first order correction in bound state per- 
turbation theory, we find all powers of a, and these must be grouped order-by-order 
with terms that appear at higher orders of bound state perturbation theory. 

The leading correction to the binding energy is O (a 4 ), and producing these correc- 
tions is a much more serious test of the renormalization procedure than the calculation 
shown above. To what order in the coupling must the Hamiltonian be computed to 
correctly reproduce all masses to (9(a 4 )? The leading error can be found in the electron 
mass itself. With the Hamiltonian given above, two-loop effects would show errors in 
the electron mass that are 0(a 2 A 2 ). This would appear to present a problem for the 
calculation of the binding energy to 0(a 2 ), but remembering that the cutoff must be 
lowered so that am 2 > A 2 > a 2 m 2 , we see that the error in the electron mass is actu- 
ally of 0(a 3+s m 2 ). This means that to compute masses correctly to 0(a A ) we would 
have to compute the Hamiltonian to 0(a 4 ), which requires a fourth-order similarity 
calculation for QED. Such a calculation has not yet been completed. However, if we 
compute the splitting between bound state levels instead, errors in the electron mass 
cancel and we find that the Hamiltonian computed to 0(a 2 ) is sufficient. 

In Ref. 10 we have shown that the fine structure of positronium is correctly 
reproduced when the first- and second-order corrections from bound state perturbation 
theory are added. This is a formidable calculation, because the exact Coulomb bound 
and scattering states appear in second-order bound state perturbation theory^] 

A complete calculation of the Lamb shift in hydrogen would also require a fourth- 
order similarity calculation of the Hamiltonian; however, the dominant contribution to 
the Lamb shift that was first computed by Bethe 36 can be computed using a Hamilto- 
nian determined to 0(a). 11 In this calculation a Bloch transformation was used rather 
than a similarity transformation because the Bloch transformation is simpler and small 
energy denominator problems can be avoided in analytic QED calculations. 

The primary obstacle to using our light-front strategy for precision QED calcula- 
tions is algebraic complexity. We have successfully used QED as a testing ground for 

^There is a trick which allows this calculation to be performed using only first-order bound state 
perturbation theory. 35 The trick basically involves using a Melosh rotation. 
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this strategy, but these calculations can be done much more conveniently using other 
methods. The theory for which we believe our methods are best suited is QCD. 

LIGHT-FRONT QCD 

This section relies heavily on the discussion of positronium, because we only require 
the QCD Hamiltonian determined to 0(a) to discuss a simple confinement mechanism 
which appears naturally in light-front QCD and to complete reasonable zeroth order 
calculations for heavy quark bound states. To this order the QCD Hamiltonian in the 
quark- ant iquark sector is almost identical to the QED Hamiltonian in the electron- 
positron sector. Of course the QCD Hamiltonian differs significantly from the QED 
Hamiltonian in other sectors, and this is essential for justifying my choice of Hq for 
non-perturbative calculations. 

The basic strategy for doing a sequence of (hopefully) increasingly accurate QCD 
bound state calculations is almost identical to the strategy for doing QED calculations. 
I use coupling coherence to find an expansion for H k in powers of the QCD coupling 
constant to a finite order. I then divide the Hamiltonian into a non-perturbative part, 
Hq , and a perturbative part, V A . The division is based on the physical argument that 
adding a parton in an intermediate state should require more energy than indicated 
by the free Hamiltonian, and that as a result these states will 'freeze out' as the cutoff 
approaches Aq C d- When this happens the evolution of the Hamiltonian as the cutoff 
is lowered further changes qualitatively, and operators that were consistently canceled 
over an infinite number of scales also freeze, so that their effects in the few parton 
sectors can be studied directly. A one-body operator and a two-body operator arise in 
this fashion, and serve to confine both quarks and gluons. 

The simple confinement mechanism I outline is certainly not the final story, but it 
may be the seed for the full confinement mechanism. One of the most serious problems 
we face when looking for non-perturbative effects such as confinement is that the search 
itself depends on the effect. A candidate mechanism must be found and then shown to 
self-consistently produce itself as the cutoff is lowered towards Aq C d- 

Once we find a candidate confinement mechanism, it is possible to study heavy 
quark bound states with little modification of the QED strategy. Of course the results 
in QCD will differ from those in QED because of the new choice of TLq , and in higher 
orders because of the gluon interactions. When we move on to light quark bound states, 
it becomes essential to introduce a mechanism for chiral symmetry breaking. 7 ' 1 I will 
discuss this briefly at the end of this section. 

When we compute the QCD Hamiltonian to O(a), several significant new features 
appear. First are the familiar gluon interactions. In addition to the many gluon interac- 
tions found in the canonical Hamiltonian, there are modifications to the instantaneous 
gluon exchange interactions, just as there were modifications to the electron-positron 
interaction. For example, a Coulomb interaction will automatically arise at short dis- 
tances. In addition the gluon self-energy differs drastically from the photon self-energy. 

The photon develops a self-energy because it mixes with electron-positron pairs, 
and this self energy is 0(aA 2 / P + ). When the cutoff is lowered below 4m 2 , this mass 
term vanishes because it is no longer possible to produce electron-positron pairs. For 
all cutoffs the small bare photon self-energy is exactly canceled by mixing with pairs 
below the cutoff. I will not go through the calculation, but because the gluon also 
mixes with gluon pairs in QCD, the gluon self-energy acquires an infrared divergence, 
just as the electron did in QED. In QCD both the quark and gluon self-energies are 
proportional to aA 2 ln(l/e)/P + , where e is the secondary cutoff on parton longitudinal 
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momenta introduced in the last section. This means that even when the primary cutoff 
A 2 is finite, the energy of a single quark or a single gluon is infinite, because we are 
supposed to let e — * 0. One can easily argue that this result is meaningless, because 
the relevant matrix elements of the Hamiltonian are not even gauge invariant; however, 
since we must live with a variational principle when doing Hamiltonian calculations, 
this result may be useful. 

In QED I argued that the bare electron self-energy was a complete lie, because 
the bare electron mixes with photons carrying arbitrarily small longitudinal momenta 
to cancel this bare self-energy and produce a finite mass physical electron. However, in 
QCD there is no reason to believe that this perturbative mixing continues to arbitrarily 
small cutoffs. There are no massless gluons in the world. In this case, the free QCD 
Hamiltonian is a complete lie and cannot be trusted at low energies. 

On the other hand, coupling coherence gives us no choice about the quark and 
gluon self-energies as computed in perturbation theory. These self-energies appear 
because of the behavior of the theory at extremely high energies. The question is not 
whether large self-energies appear in the Hamiltonian. The question is whether these 
self-energies are canceled by mixing with low energy multi-gluon states. I argue that 
this cancellation does not occur, and that the infrared divergent quark and gluon self- 
energies should be included in Hq . The transverse scale for these energies is the running 
scale A, and over many orders of magnitude we should see the self-energies canceled by 
mixing. However, as the cutoff approaches Aqcd-, I speculate that these cancellations 
cease to occur because perturbation theory breaks down and a mass gap between states 
with and without extra gluons appears. 

But if the quark and gluon self-energies diverge, and the divergences cannot be 
canceled by mixing between sectors with an increasingly large number of partons, how is 
it possible to obtain finite mass hadrons? The parton-parton interaction also diverges, 
and the infrared divergence in the two-body interaction exactly cancels the infrared 
divergence in the one-body operator for color singlet states. 

Of course, the cancellation of infrared divergences is not enough to obtain con- 
finement. The cancellation is exact regardless of the relative motion of the partons 
in a color singlet state, and confinement requires a residual interaction. I will show 
that the 0(a) QCD Hamiltonian produces a logarithmic potential in both longitudinal 
and transverse directions. I will not discuss whether a logarithmic confining poten- 
tial is 'correct,' but to the best of my knowledge there is no rigorous demonstration 
that the confining interaction is linear, and a logarithmic potential may certainly be 
of interest phenomenologically for heavy quark bound states. 37 ' 38 I would certainly be 
delighted if a better light-front calculation produces a linear potential, but this may 
not be necessary even for light hadron calculations. 

The calculation of how the quark self-energy changes when a similarity transforma- 
tion lowers the cutoff on energy transfer is almost identical to the electron self-energy 
calculation. Following the steps in the section on positronium, we find the one-body 
operator required by coupling coherence, 

e£*(p) 



where C F = (N 2 - 1)/(2N) for a SU(N) gauge theory. 

The calculation of the quark-antiquark interaction required by coupling coherence 
is also nearly identical to the QED calculation. Keeping only the infrared singular parts 
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of the interaction, as was done in QED, 
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The instantaneous gluon exchange interaction is 

vLtam = -^glc F \jptptkt H 5 aia2 5 Xl X 2 

x e(\q+\ - eP + ) e(^- | V - x + fcr -pj - fc 2 " |^ . (155) 

Just as in QED the coupling coherent interaction induced by gluon exchange above the 
cutoff partially cancels instantaneous gluon exchange. For the discussion of confinement 
the part of Vcoh that remains is not important, because it produces the short range part 
of the Coulomb interaction. However, the part of the instantaneous interaction that is 
not canceled is 
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Note that this interaction contains a cutoff that projects onto exchange energies 
below the cutoff, because the interaction has been screened by gluon exchange above the 
cutoffs. This interaction can become important at long distances, if parton exchange 
below the cutoff is dynamically suppressed. In QED I argued that this singular long 
range interaction is exactly canceled by photon exchange below the cutoff, because 
such exchange is not suppressed no matter how low the cutoff becomes. Photons are 
massless and experience no significant interactions, so they are exchanged to arbitrarily 
low energies as effectively free photons. This cannot be the case for gluons. 

For the discussion of confinement, I will place only the most singular parts of the 
quark self-energy and the quark- ant iquark interaction in TCq. To see that all infrared 
divergences cancel and that the residual long range interaction is logarithmic, we can 
study the matrix element of these operators for a quark-antiquark state, 
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where r and s are color indices and I will choose to be a color singlet and drop color 
indices. The cancellations we find do not occur for the color octet configuration. The 
matrix element is, 

(${P)\H%\V(P')) = 16tt 3 P+5 3 (P - P') x 

dyd 2 K 
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Here I have chosen a frame in which the center-of-mass transverse momentum is zero, 
assumed that the longitudinal momentum scale introduced by the cutoffs is that of the 
bound state, and used Jacobi coordinates, 
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The first thing I want to do is show that the last term is divergent and the di- 
vergence exactly cancels the first term. My demonstration is not elegant, but it is 
straightforward. The divergence results from the region y\ ~ y-i- In this region the 
second and third cutoffs restrict {k\ — k 2 ) 2 to be small compared to A 2 , so we should 
change variables, 
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Using these variables we can approximate the above interaction near q = and y — 0. 
The double integral becomes 
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where i] is an arbitrary constant that restricts \y\ from becoming large. Completing the 
q and y integration we get 
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The divergent part of this exactly cancels the first term on the right-hand side 
of Eq. (122). This cancellation occurs for any state, and this cancellation is unusual 
because it is between the expectation value of a one-body operator and the expecta- 
tion value of a two-body operator. The cancellation resembles what happens in the 
Schwinger model and is easily understood. It results from the fact that a color singlet 
has no color monopole moment. If the state is a color octet the divergences are both 
positive and cannot cancel. Since the cancellation occurs in the matrix element, we can 
let e — > before diagonalizing TCq . 

The fact that the divergences cancel exactly does not indicate that confinement 
occurs. This requires the residual interactions to diverge at large distances, which 
means small momentum transfer. Equivalently, we need the color dipole self-energy 
to diverge if the color dipole moment diverges because the partons separate to large 
distance. 

My analysis of the residual interaction is neither elegant nor complete. For a more 
complete analysis see Ref.12. I show that the interaction is logarithmic in the longitudi- 
nal direction at zero transverse separation and logarithmic in the transverse direction at 
zero longitudinal separation, and I present the full angle-averaged interaction without 
derivation. 

In order to avoid the infrared divergence, which is canceled, I compute spatial 
derivatives of the potential. First consider the potential in the longitudinal direction. 
Given a momentum space expression, set x± = and the Fourier transform of the 
longitudinal interaction requires the transverse momentum integral of the potential, 



§-V(z) = J ^L iqz V(q ± , q z )e^ z . (163) 



We are interested only in the long range potential, so we can assume that q z is arbitrarily 
small during the analysis and approximate the step functions accordingly. 
For our interaction this leads to 
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Completing the q±_ integration we have 
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To see that the term involving a cosine in the next-to-last line produces a short range 
potential, simply integrate it. At large |:r~|, which is the only place we can trust 
our approximations for the original integrand, this yields a logarithmic potential, as 
promised. 
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Next consider the potential in the transverse direction. Here we can set x = 
and get 

d . 2 , f dq + d 2 q 



dx ± 
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Here I have used the fact that the integration is dominated by small q + to simplify the 
integrand again. Completing the q + integration this becomes 
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7T 2 x\ 
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Once again, this is the derivative of a logarithmic potential, as promised. The strength 
of the long-range logarithmic potential is not spherically symmetrical in these coor- 
dinates, with the potential being larger in the transverse than in the longitudinal di- 
rection. Of course, there is no reason to demand that the potential is rotationally 
symmetric in these coordinates. 

The angle-averaged potential for two heavy quarks with mass M is 12 

g 2 A C F A 2 / Si(n) (l-cos(ft)) sin(ft) 5 x , . 

(ln(tt) - Ci(fc) + 2-^ - * ^A-^ + — ^ - - + lE ) , (168) 

where, 

A 2 
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I have assumed the longitudinal momentum scale in the cutoff equals the longitudinal 
momentum of the state. 

The full potential is not naively rotationally invariant. Of course, the two body 
interaction in QED is also not rotationally invariant. The leading term in an expansion 
in powers of momenta yields the rotationally invariant Coulomb interaction, but higher 
order terms are not rotationally invariant. These higher order terms do not spoil 
rotational invariance of physical results because they must be combined with higher 
order interactions in the effective Hamiltonian, and with corrections involving photons 
at higher orders in bound state perturbation theory. There is no reason to expect, and 
no real need for a potential that is rotationally invariant; however, to proceed we need 
to decide which part of the potential must be treated non-perturbatively. The simplest 
reasonable choice is the angle-averaged potential, which is what we have used in heavy 
quark calculations. 12 ' 13 

Had we computed the quark-gluon or gluon-gluon interaction, we would find es- 
sentially the same residual long range two-body interaction in every Fock space sector, 
although the strengths would differ because different color operators appear. In QCD 
gluons have a divergent self-energy and experience divergent long range interactions 
with other partons if we use coupling coherence. In this sense, the assumption that 
gluon exchange below some cutoff is suppressed is consistent with the Hamiltonian 
that results from this assumption. To show that gluon exchange is suppressed when 
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A — > Aqcd, rather than some other scale (i.e., zero as in QED), a non-perturbative 
calculation of gluon exchange is required. This is exactly the calculation bound state 
perturbation theory produces, and bound state perturbation theory suggests how the 
perturbative renormalization group calculation might be modified to generate these 
confining interactions self-consistently. 

If perturbation theory, which produced this potential, continues to be reasonable, 
the long range potential will be exactly canceled in QCD just as it is in QED. We need 
this exact cancellation of new forces to occur at short distances and turn off at long dis- 
tances, if we want asymptotic freedom to give way to a simple constituent confinement 
mechanism. At short distances the divergent self-energies and two-body interactions 
cannot be ignored, but they should exactly cancel pairwise if these divergences appear 
only when gluons are emitted and absorbed in immediately successive vertices, as I 
have speculated. The residual interaction must be analyzed more carefully at short 
distances, but in any case a logarithmic potential is less singular than a Coulomb in- 
teraction, which does not disturb asymptotic freedom. This is the easy part of the 
problem. The hard part is discovering how the potential can survive at any scale. 

A perturbative renormalization group will not solve this problem. The key I have 
suggested is that interactions between quarks and gluons in the intermediate states 
required for cancellation of the potential will eventually produce a non-negligible energy 
gap. I am unable to detail this mechanism without an explicit calculation, but let me 
sketch a naive picture of how this might happen. 

Focus on the quark-antiquark-gluon intermediate state, which mixes with the 
quark- ant iquark state to screen the long range potential. The free energy of this in- 
termediate state is always higher than that of the quark-antiquark free energy, as is 
shown using a simple kinematic argument. 1 However, if the gluon is massless, the en- 
ergy gap between such states can be made arbitrarily small. As we use a similarity 
transformation to run a vertex cutoff on energy transfer, mixing persists to arbitrarily 
small cutoffs since the gap can be made arbitrarily small. Wilson has suggested using a 
gluon mass to produce a gap that will prevent this mixing from persisting as the cutoff 
approaches 1 Aqcd, and I am suggesting a slightly different mechanism. 

If we allow two-body interactions to act in both sectors to all orders, even the 
Coulomb interaction can produce quark-antiquark and quark-antiquark-gluon bound 
states. In this respect QCD again differs qualitatively from QED because the photon 
cannot be bound and the energy gap is always arbitrarily small even when the electron- 
positron interaction produces bound states. If we assume that a fixed energy gap arises 
between the quark-antiquark bound states and the quark-antiquark-gluon bound states, 
and that this gap establishes the important scale for non-perturbative QCD, these states 
must cease to mix as the cutoff goes below the gap. 

An important ingredient for any calculation that tries to establish confinement self- 
consistently is a seed mechanism, because it is possible that it is the confining interaction 
itself which alters the evolution of the Hamiltonian so that confinement can arise. I 
have proposed a simple seed mechanism whose perturbative origin is appealing because 
this allows the non-perturbative evolution induced by confinement to be matched on 
to the perturbative evolution required by asymptotic freedom. 

I provide only a brief summary of our heavy quark bound state calculations, 12, 13 
and refer the reader to the original articles for details. We follow the strategy that 
has been successfully applied to QED, with modifications suggested by the fact that 
gluons experience a confining interaction. We keep only the angle-averaged two-body 
interaction in T~Lq , so the zeroth order calculation only requires the Hamiltonian matrix 
elements in the quark-antiquark sector to 0(a). All of the matrix elements are given 
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above. 

For heavy quark bound states 13 we can simplify the Hamiltonian by making a 
nonrelativistic reduction and solving a Schrodinger equation, using the potential in Eq. 
(168). We must then choose values for A, a, and M. These should be chosen differently 
for bottomonium and charmonium. The cutoff for which the constituent approximation 
works well depends on the constituent mass, as in QED where it is obviously different 
for positronium and muonium. In order to fit the ground state and first two excited 
states of charmonium, we use A = 2.5GeV, a = 0.53, M c = 1.6GeV. In order to 
fit these states in bottomonium we use A = A.9GeV, a = 0.4, and Mb = A.8GeV.f\ 
Violations of rotational invariance from the remaining parts of the potential are only 
about 10%, and we expect corrections from higher Fock state components to be at least 
of this magnitude for the couplings we use. 

These calculations show that the approach is reasonable, but they are not yet very 
convincing. There are a host of additional calculations that must be done before the 
success of this approach can be judged, and most of them await the solution of several 
technical problems. 

In order to test the constituent approximation and see that full rotational invari- 
ance is emerging as higher Fock states enter the calculation, we must be able to include 
gluons. If gluons are massless (which need not be true since the gluon mass is a relevant 
operator that must in principle be tuned), we cannot continue to employ a nonrelativis- 
tic reduction. In any case, we are primarily interested in light hadrons and must learn 
how to perform relativistic calculations to study them. The primary difficulty is that 
evaluation of matrix elements of the interactions given above involves high dimensional 
integrals which display little symmetry in the presence of cutoffs. Worse than this is 
the fact that the confinement mechanism requires cancellation of infrared divergences 
which have prevented us from using Monte Carlo methods to date. 31 These difficulties 
are avoided when a nonrelativistic reduction is made, but there is little more that we 
can do for which such a reduction is valid. 

I conclude this section with a few comments on chiral symmetry breaking. While 
quark-quark, quark-gluon, and gluon-gluon confining interactions appear in the hamil- 
tonian at 0(a), chiral symmetry is not broken at any finite order in the perturbative 
expansion; so symmetry-breaking operators must be allowed to appear ab initio. Light- 
front chiral symmetry differs from equal-time chiral symmetry in several interesting and 
important aspects. 1 ' 7> 9 For example, quark masses in the kinetic energy operator do not 
violate light-front chiral symmetry; and the only operator in the canonical hamiltonian 
that violates this symmetry is a quark-gluon coupling linear in the quark mass. Again, 
the primary technical difficulty is the need for relativistic bound state calculations, 
and real progress on chiral symmetry cannot be made before we are able to perform 
relativistic calculations with constituent gluons. If transverse locality is maintained, 
simple renormalization group arguments indicate that chiral symmetry should be bro- 
ken by a relevant operator in order to preserve well-established perturbative results at 
high energy. The only relevant operator that breaks chiral symmetry involves gluon 
emission and absorption, leading to the conclusion that the pion contains significant 
constituent glue. This situation is much simpler than what we originally envisioned, 1 
because it does not require the addition of any operators that cannot be found in the 
canonical QCD hamiltonian; and we have long known that relevant operators must be 
tuned to restore symmetries. 



"There are several minor errors in Ref. 13, which are discussed in Ref.14. I have also chosen V + = P + , 
and in those papers V + = P + /2; and I have absorbed this change in a redefinition of A. 
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